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Abstract 

An  experimental  verification  of  the  end  effect  theory  of  long  cylindrical 
Langmuir  probes  in  a high-speed  collisionless  plasma  flow  is  reported.  This  effect 
is  observed  as  a sharp  departure  from  the  predicted  'infinite  probe'  current  when 
a strongly  negative  probe  is  aligned  with  the  flow.  Both  the  height  and  the  width 
of  the  ion  current  peak  measured  about  the  zero-degree  attitude  are  strongly 
sensitive  to  ion  temperature . The  effect  is  verified  for  a mesothermal  argon 
plasma  with  probe  radius  to  Debye  length  ratio  0.005  < e < 0.5  and  for  ion  to 
electron  temperature  ratios  close  to  unity.  A scheme  is  developed  to  obtain  a 
self-consistent  measurement  of  the  ion  temperature  using  the  basic  characteristics 
of  the  finite-probe  angular  response  after  a modification  of  J.  R.  Sanmartin's 
theory. 

A versatile  ion  source  was  constructed  after  a state-of-the-art  electro- 
static ion  thruster.  It  provides,  together  with  a large  size  vacuum  chamber,  a 
simulation  facility  for  which  the  plasma  stream  characteristics  can  be  tailored 
and  the  ion  speed  ratio  easily  modeled  for  diagnostics  purposes. 

A range  of  flow  conditions  and  of  probe  parameters  are  investigated 
and  the  overall  applicability  and  reliability  of  the  method  to  measure  T^  are 
assessed.  Potential  areas  of  applications  are  outlined. 
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Sommaire 


Un  effet  de  bout  de  longues  sondes  cylindriques  de  Langmuir  a 4t4 
v4rifi4  exp^rimentalement  pour  un  flux  molbculaire  de  plasma  & grande  vitesse 
d'^coulement.  L' effet  observe  se  rkvfele  trbs  fortement  en  des accord  avec  la 
thborie  bien  connue  des  sondes  cylindriques  "infinies",  quand  la  sonde  est 
par  alible  au  faisceau.  La  hauteur  et  la  dimension  angulaire  du  pic  de  courant 
ionique  observe,  sont  toutes  deux  trks  sensibles  k la  temperature  des  ions. 
L'effet  de  bout  a 4t4  v4rifi4  k l'aide  d'un  plasma  d'argon  mesothermal,  ok 
les  temperatures  ionique  et  4lectronique  4taient  k peu  prks  dgales,  et  pour 
des  rapports  rayon  de  sonde- longueur  de  Debye,  denotes  e,  allant  de  0.005  k 
0.5  environ.  Une  m4thode  autocritique  de  mesure  de  la  temperature  ionique, 
k partir  des  donn4es  caract4ristiques  de  la  response  courant-angle  d'attaque, 
a 4t4  ddmontrbe  en  utilisant  la  thborie  quelque  peu  modifi4e  de  J.  R.  Sanmartin. 
Une  source  d'ions  d1 operation  trks  flexible  inspir4e  d'un  prototype  de  propul- 
seur  ionique  a 4t4  construite  pour  cette  etude.  L'accbl^rateur  d'ions  ainsi 
qu’un  caisson-k-vide  de  grandes  dimensions  constituent  une  enceinte  de 
simulation  d'ecoulements  de  plasma  trks  flexible.  Les  caracteristiques 
principales  ainsi  que  les  proprikt^s  de  vitesse  d'ecoulement  et  de  vitesse 
thermique  de  ceux-ci  peuvent  ktre  decrites  k l'aide  d'un  modkle  single,  et 
modifiees  k des  fins  experimentales . Toute  une  gamme  de  conditions  d'4coule- 
ment  et  de  caracteristiques  de  sonde  ont  4t4  4tudi4es.  L'utilisation  ainsi 
que  la  precision  de  la  mkthode  de  mesure  de  la  temperature  ionique  ont  4t4 
dkfinies  et  quelques  applications  y sont  mentionn^es . 
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Introductory  Comments 


The  study  of  the  behaviour  of  electrically  biased  metallic  probes  in 
ionized  media  has  a considerable  bearing  on  the  simulation  and  on  the  diagnostics 
of  space-flight  environment  and  of  laboratory  plasmas  of  many  kinds.  Since  the 
original  work  of  Langmuir,  Mott-Smith  and  Tonks  in  the  mid- twenties  (Refs.  1,  2) 
an  inpressive  number  of  contributions  to  the  literature  have  appeared  on  this 
subject.  A prime  motivation  is  invariably  to  find  new  analytical  solutions  to 
the  fundamental  problems  associated  with  the  collection  of  charged  particles  under 
a growing  variety  of  conditions.  In  particular,  the  problems  concerning  flowing 
plasmas  of  many  origins  have  led  to  the  refining  of  numerous  schemes  for  extracting 
more  and  more  parameters  directly  from  the  characteristics  of  the  current  collec- 
tion characteristics  - namely,  ion  and  electron  density,  electron  temperature, 
collisional  parameters,  chemical  and  thermodynamic  properties,  particle  distribu- 
tion function  and  generally  speaking  with  greater  difficulty,  ion  temperature. 

These  efforts  have  been  devoted  to  the  study  of  plasmas  of  natural  sources  (iono- 
sphere, sun  and  stars,  interplanetary  space,  lightning,  etc.)  and  now  more  than 
ever  to  plasmas  of  man-made  origins  such  as  those  of  controlled  thermonuclear 
fusion  machines,  of  communication  devices,  those  of  electromagnetic  propulsion 
devices,  of  laser-matter  interactions  and  of  plasma  generators  for  large  energy 
c onver  si on  sy s t ems . 

In  recent  years  the  work  of  Laframboise  (Ref.  3)  who  developed  an 
analysis  based  on  the  related  work  of  Bernstein  and  Rabinowitz  (Ref.  4),  has 
been  recognized  and  established  as  a standard  for  two-component  collisionless 
Maxwellian  plasmas  at  rest;  it  constitutes  in  fact  a firm  basis  to  draw  from  in 
studies  of  dynamic  and  unsteady  cases. 

It  has  been  the  contention  of  several  users  of  electrostatic  probes, 
often  simply  called  Langmuir  probes,  in  flowing  collisionless  plasma  streams 
that  the  flow  does  not  significantly  alter,  if  at  all,  the  response  of  suf ficiently 
long  cylindrical  probes  whether  the  probe  axis  is  oriented  parallel  or  perpendicular 
to  the  flow  direction  (Refs.  5-13)*  On  the  basis  of  such  arguments  it  is  then 
plausible  to  obtain  a successful  interpretation  of  the  probe  response  in  terms  of 
the  exhaustive  theories  available  for  'static'  cases.  Therefore  a local  observation 
of  the  main  plasma  parameters  is  available  without  the  need  for  a special  theory 
for  flowing  conditions.  These  premises  have  in  fact  enabled  numerous  workers  to 
verify  the  'static  plasma'  approach  in  supersonic  flows  with  large  probe  radius 
to  Debye  length  ratios  (Refs.  l4-l6) . 

However,  early  experimental  results  reported  by  L.  H.  Brace  (Ref.  l6a, 
1965)  conducted  with  long  cylindrical  probes  on-board  the  Explorer-17  satellite, 
have  shown  a markedly  large  attitude  sensitivity  of  the  current  collection  when 
the  probes  were  nearly  aligned  with  the  satellite  velocity  vector.  In  fact  a 
sharp  ion  current  peak  was  observed  about  the  0°  attitude  whereas  the  so-called 
'orbital  limited'  collection  expected  would  have  shown  the  slow  monotonic  decrease 
which  characterizes  the  'infinite'  probe  behaviour.  Bettinger  and  Chen  were  the 
first  to  analyze  these  findings  (Ref.  17).  They  explained  Brace's  observations 
in  terms  of  an  end  effect  of  long  cylindrical  probes  when  the  Debye  length  is 
large  compared  with  the  probe  radius  and  when  the  flow  speed  is  high,  even  though 
probe  aspect  ratios  ($/ d)  of  nearly  400  were  considered. 

This  anomalous  behaviour  obviously  in  contradiction  with  the  'infinite' 
probe  expectations  was  subsequently  observed  in  a laboratory  simulation  of  a 


similar  ionospheric  environment  by  Hester  and.  Sonin  in  1969  (Refs.  18,20).  They 
have  published  a qualitative  assessment  of  this  end  effect  based  on  a numerical 
study  in  order  to  reconcile  the  theory  with  their  own  observations;  they  also 
indicated  a fundamental  disagreement  with  Bettinger  and  Chen's  explanations. 

Recently  the  work  of  Sanmartin  (1971)  has  finalized  the  assessment  of  the  effect 
in  an  extensive  analysis  (Refs.  21-23).  As  it  was  originally  pointed  out  by 
Hester  and  Sonin  the  finite-probe  theory  unveiled  a considerable  interest  since 
it  affected  the  entire  conceptual  interpretation  of  the  probe  response  in  flowing 
plasmas,  yet  more  so  because  the  effect  appeared  to  be  sensitive  to  some  plasma 
parameters  under  a wide  range  of  stream  conditions.  In  fact  they  alluded  to 
the  possibility  of  using  the  effect  for  diagnostic  purposes  - notably,  to  extract 
some  valuable  information  from  the  strong  attitude  sensitivity  and  possibly  from 
the  ion  temperature  sensitivity  contained  in  the  very  structure  of  current  angular 
response.  This  last  point  was  reiterated  by  Sanmartin  (1972)  who  elaborated  the 
foundations  of  a simple  method  to  obtain  the  ion  temperature  from  the  ion  peaks 
observed.  Shortly  after,  Fournier  and  Pigache  (Ref.  57)  have  published  some 
quantitative  attests  to  apply  Sanmartin' s analytical  description.  More  specifically 
they  verified  the  end  effect  theory  using  a long  cylindrical  probe  within  a plasma 
stream,  under  simulated  ionospheric  conditions,  and  they  also  reported  some  estimates 
of  the  ion  speed  ratio. 

The  many  aspects  of  the  end  effect  of  cylindrical  probes  have  generated 
a considerable  interest  in  the  UTIAS  Low  Density  Plasma  Dynamics  Group,  at  the 
onset  on  this  research  program,  especially  so  because  of  the  sharp  attitude 
sensitivity  it  entailed.  At  the  time,  in  early  1972,  an  experimental  investigation 
was  underway  to  study  the  exhaust  characteristics  of  a small-scale  electrostatic 
ion  thruster.  A bench  prototype  which  was  built  and  tested  by  the  author  in 
this  laboratory  was  used  to  reproduce  the  beam  characteristics  of  state-of-the-art 
ion  thrusters,  such  as  the  ones  proposed  for  station-keeping  and  attitude- control 
of  synchronous  satellites.  The  need  for  a very  fine  pointing  accuracy  was 
desirable,  in  particular  for  the  study  of  the  beam  vectoring  ability  of  flight- 
worthy propulsion  units. 

Actually  the  use  of  small  ion  accelerators  of  that  type  in  large 
simulation  vacuum  chambers  lends  itself  quite  naturally  to  a number  of  fundamental 
studies  involving  high-speed  plasma  flows,  such  as  the  modelling  of  ionospheric 
sounding  conditions  for  instance  (Refs.  20,29,30).  Moreover,  since  the  end  effect 
also  showed  the  promising  advantages  of  a simple  diagnostic  tool  to  measure  the 
ion  temperature  in  fast  drifting  plasmas  we  shifted  our  prime  motivation  and  we 
set  out  to  establish  the  basis  of  a systematic  study  of  the  anomalous  probe 
response.  The  main  research  objectives  which  were  elaborated  in  this  respect, 
and  which  find  their  conclusion  in  the  present  report  are  summarized  as  follows: 

1.  Provide  a facility  to  achieve  collisionless-high-speed  plasma  con- 
ditions suitable  for  laboratory  simulation  purposes  and  for  a variety  of  plasma 
investigations.  The  stream  parameters  should  be  controllable,  easily  monitored 
and  systematically  tailored  to  satisfy  a wide  range  of  experimental  requirements 
and  to  allow  accurate  beam  diagnostics  (see  Chapter  n) . 

2.  Verify  the  salient  features  of  the  end  effect  of  cylindrical  probes 
described  in  the  theory  proposed  by  J.  R.  Sanmartin,  for  a wide  range  of  flow 
conditions  and  of  probe  characteristics.  It  is  also  necessary  to  describe  and 

to  validate  the  analysis  of  the  effect  using  a set  of  nondimensional  similarity 
parameters  independent  of  electron  temperature  as  opposed  to  the  one  originally 
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selected  by  Hester  and  Sonin.  In  order  to  conplete  the  experimental  results 
provided  by  these  authors  and  later  by  Fournier  and  Pigache,  a systematic 
investigation  is  required.  In  particular  the  reliability  of  the  theory  should 
be  checked  (Chapter  I)  with  a drifting  plasma  in  which  the  ion  temperature  is  of 
the  same  order  as  the  electron  temperature,  since  only  the  limiting  cold-ion 
case,  T^/Te  ~ 0,  has  been  implemented.  Moreover  the  large-sheath  regime  should 
be  studied  down  to  values  of  probe  radius  approximately  equal  to  electron  Debye 
length  (rp/-Aj)  ~ l)  in  order  to  cover  both  weak  and  strong  probe  effects. 

3.  Develop  a method,  using  Sanmar tin's  analysis  which  would  demonstrate 
the  ability  of  the  end  effect  to  serve  as  a diagnostic  tool  in  obtaining  the  ion 
temperature  by  direct  measurement  from  the  probe  angular  response  (see  Chapter  IV). 
More  specifically,  a flow-field  model  would  be  desirable  to  provide  the  necessary 
grounds  for  collaring  the  theoretical  predictions  with  the  estimates  of  the  ion 
speed  ratio  and  the  ion  temperature  derived  by  direct  measurement  (see  Chapter 
III) . 

4.  Draw  the  appropriate  conclusions  concerning  the  range  of  validity 

of  the  theory  and  the  reliability  of  the  method  to  measure  Tj_  in  order  to  delineate 
some  possible  areas  of  applications. 
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I.  THEORETICAL  CONSIDERATIONS 


1.1  Pertinent  Theory  on  the  End  Effect 

1.1.1  Statement  of  the  Problem 


The  problem  at  hand  is  that  of  the  collection  of  positive  ions  by  a 
cylindrical  probe  in  a high-speed,  steady  and  collisionless  plasma  flow,  for  the 
regime  where  the  sheath  radius  is  large  with  respect  to  the  probe  radius.  At  the 
low  densities  of  interest  (A  » r)  the  sheath  radius  is  basically  of  the  same 
size  as  the  electron  Debye  length  the  regime  which  is  appropriate  to  our 
problem  is  the  one  described  by  the  hatched  region  in  Fig.  1.1.  For  simplicity 
we  consider  a plasma  composed  of  a single  species  of  singly- charged  positive  ions 
of  mass  M,  at  temperature  T^,  and  of  electrons  at  temperature  Te.  The  flow 
conditions  are  those  of  a drifting  mesothermal  plasma  (ne  = n^),  i.e. 


(1.1) 


and  the  probe  bias  0p,  sufficiently  negative  to  repel  effectively  all  electrons, 
is  large  enough  so  that  electron  density  distribution  is  accurately  described 
by  a Boltzmann  distribution  and 


kTe  « |e0p  | (1.2) 

The  probe  geometry  is  depicted  in  Fig.  1.2  where  i is  the  probe  length 
and  d its  diameter.  We  have  the  probe  condition 


( 


» 1 


where  the  sheath  radius  (a)  is  basically  given  by 


(1.3) 


a 


Under  these  conditions  the  end  effect  of  a cylindrical  Langmuir  probe  is  observed 
as  a sharp  increase  of  the  ion  current  collected,  relative  to  the  ion  current 
expected  using  the  'infinite'  probe  theory,  when  the  probe  is  aligned  with  the 
flow.  The  current  predicted  for  the  'infinitely'  long  cylinder  under  similar 
conditions  is  simply  the  so-called  'orbital-limited'  current  derived  by  Langmuir 
and  Mott- Smith  (Ref.  l)  and  which  is  given  by 

J»  = en^  u • d*i  • [ sin20  + <fcp]5  (1.4) 

where  in  this  report  we  choose  to  use  a nondimen si onal  probe  bias  <$p  given  by 
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e(0  - <P  ) 

<*  ^ __2 SL 

^ | Mu2 


(I. 


i.e.,  we  use  the  kinetic  energy  of  the  ion  motion  to  nondimension&liz*  ■ r > 


the  more  conventional  electron  thermal  energy  kTe. 
at  the  parallel  probe  (0  s 0)  is 


Therefore  the  curr-ii.t  le: 


' bAf 


en,  u 


(9  - 0) 


-T-^v 


The  ion  current  to  a probe  of  finite  length  in  fact  closely  t - 
'infinite'  probe  description  only  at  high  angles  of  incidence.  As  the  j i -t-  » 

brought  into  alignment  with  the  flow  a striking  deviation  occurs  in  th-  f t. 
sharp  peak  in  the  ion  current,  as  shown  in  Fig.  1.3  for  four  arbitrary  ntfeie.w 
experimental  conditions. 


This  current  peak  results  from  the  finite  length  of  the  probe, 
is  considered  as  an  end  effect.  The  probe  is  essentially  infinite  in’  t-rnu  : 
flow  condition  if  its  nondimensional  length,  defined  by  t» 


(1.7) 


according  to  the  work  of  Hester  and  Sonin  (Ref.  18) , approaches  infinity.  Briefly, 
the  physical  explanation  of  the  peak  is  that  it  represents  the  current  from  those 
ions  which  enter  the  field  of  the  probe  through  the  front  end  of  the  cylindrical 
sheath  which  surrounds  it,  rather  than  through  the  sides  as  for  an  'infinitely' 
long  cylinder.  With  the  sheath  radius  large  compared  with  the  probe  radius  and  the 
bulk  speed  large  compared  with  the  ion  random  speed,  the  ion  flux  through  the  front 
of  the  sheath  can  be  significant  when  the  probe  is  aligned  with  the  flow.  However, 
when  the  probe  is  turned  by  a small  angle  virtually  all  of  the  ions  entering  the 
front  'scoop'  have  sufficient  angular  momentum  transverse  to  the  probe  to  escape 
the  field,  and  hence  the  sharp  drop  in  current  registered  toward  the  smaller  corre- 
sponding 'infinite'  probe  value.  It  also  can  be  seen  intuitively  how  the  end  effect 
can  be  a strong  function  of  the  ion  temperature,  since  some  thermal  energy  may  give 
a large  fraction  of  the  ion  population  entering  the  front  of  the  sheath  sufficient 
angular  momentum  to  escape  the  probe,  even  at  zero  degree  of  incidence. 


1.1.2  Analytical  Description 

The  analysis  reported  in  this  section  reproduces  the  relevant  elements 
from  the  work  of  Sanmartin  reported  in  Refs.  22  and  23.  Sanmar tin’s  approach  to 
the  problem  outlined  in  the  previous  section  is  based  on  the  following  premises: 

(i)  In  order  to  obtain  the  correct  description  of  the  ion  current 
collected  by  the  probe,  one  requires  to  know  the  exact  value 
of  the  electric  field  along  the  ion  streamlines.  However  it 
is  possible  to  derive  an  accurate  expression  for  the  electric 
field  (d0/dr)  for  a limited  radial  region  r < rm  without  having 
to  solve  simultaneously  all  three  equations  which  in  cylindrical 
coordinates  can  be  written:  (using  r,  £,  z) 


2 


where 

An  = (Znt  - ng) 
and  the  boundary  values  are 


0 (r  -»»)  =0;  0(r  = r ) = ^ 

(ii)  Once  the  electric  field  is  known  at  each  point,  the  ion 
trajectories  can  be  confuted  explicitly. 

(iii)  For  certain  conditions  relatively  fulfilled  for  the  probe 
potential,  this  steady  problem  can  be  described  with 
satisfactory  accuracy  by  an  equivalent  time- dependent 
problem  in  which  time  represents  the  flow  distance  along 
the  probe  measured  from  the  tip. 

(iv)  The  ion  current  to  the  probe  is  linear-  in  the  unperturbed 
ion  distribution  function  at  t = 0 which  is  that  of  the 
ions  just  arriving  at  the  front  end  of  the  sheath. 

The  extensive  analytical  description  which  we  only  summarize  here  is  treated 
exhaustively  in  Ref.  23.  The  corner-stone  of  Sanmartin's  arguments  as  described 
in  (i)  and  (iii)  is  justified  as  follows.  First  a list  of  reasonable  assumptions 
is  invoked,  namely,  that 

the  velocity  of  the  ions  parallel  to  the  -probe's  axis  is  basically_ 
the  unperturbed  bulk  velocity  u since  u » ^2kTj_/M  and  e0^  «|-  Mu2, 

the  front  sheath  effects  are  neglected  such  that  the  ion  distribution 
function  of  the  ions  at  t = 0 is  a drifting  Maxwellian  with  temperature 
Ti> 


3 


asymmetry  in  the  sheath  due  to  small  angles  of  attack  is  neglected 
such  that  changes  in  the  ion  angular  momentum  and  in  the  energy  can 
be  ignored, 


the  probe  is  perfectly  absorbing  for  both  species,  that  is,  the 
charged  particles  are  'annihilated'  upon  arrival  at  the  probe. 


Secondly,  upon  integrating  Poisson's  equation  obtained  from  condition  (i)  which 
simply  reads 


i | -(  T M±l'\  ~ eAn(t) 

r dr  V or  J ~ e 


we  can  obtain  the  time- dependent  electric  field.  If  one  assumes  that  An  is 
independent  of  r for  a small  region  of  interest  near  the  probe,  then 


JZ  \ M t 


t]_  + eAn(t) 

r 2e 

r=r  o 

P 


r2  - r 2 


(1.8) 


If  one  defines  the  time-dependent  field  at  the  probe  (r=r  ) by  separating  the  time 
independent  part  P 


Mill 

Sr 


(1-9) 


and  justifying  the  use  of  an  average  value  for  5(t)  and  An(t),  6 and  An 
respectively  (Ref.  23),  within  a region 


r < r < r 

p — - m 


we  get 


v*  80  -t—  / r2  - r 2 

& ~ . _R  + eAn  f £_ 

Jr  r 2e  V r 
o v 


(1.10) 


which  is  the  required  electric  field.  According  to  (ii)  therefore,  the  current 
density  j(v[_*)  can  be  computed  explicitly  and  since  it  is  linear  in  the  unper- 
turbed Maxwellian  distribution  function  (z  = -oo,  t < 0) 


f = f (v,*) 


we  have 


j(t  > 0)  mj  fo(v±*)  (J(v±*))  dVj 


(1.11) 
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where  the  * refers  to  the  simplest  form  of  the  distribution,  namely  that  for  which 
all  the  ions  have  exactly  the  same  transverse  velocity  v_j*  (magnitude  and  direction) 
along  an  arbitrary  polar  angle  £.  Thus  the  task  of  getting  j(t  > 0),  the  probe 
current  density,  is  reduced  to  that  of  evaluating 


i*  = 


i.e.  using  the  appropriate  form  for  fQ,  I. 11  becomes 


J(e)  = ^f  a cf  dv±*  ~ exp  |u±  - V,.} 

0 0 1 ci 


i* 


(1.12) 


where  £ is  the  arbitrary  angle  between  Vj*  and  the  normal  to  the  probe  and  where 
0 is  the  angle  between  the  probe  and  the  "bulk  velocity  u such  that 


Uj_  = u sine 


and 

2 

luj,  - Vj*|  = vx*2  + (u  sin0')2  - 2Vj*  u sin0  cos£ 


The  exact  analytical  representation  of  Eq.  1.12  is  by  no  means  straightforward  to 
obtain,  and  was  evaluated  by  Sanmartin  in  Ref.  22.  Using  the  ’infinite'  probe 
current  value  at  0 = 0 obtained  from  1.6  to  nondimensionalize  1.12,  Sanmartin 
gets 


f-  s T"  (B,  t e,  P,  Ma,  e)  (1.13) 

U00  ^00  * 

After  a close  inspection  of  the  conqplete  analytical  expression  for  1.13  one  can 
notice  that  the  same  result  can  be  reproduced  exactly  by  choosing  to  use  a set 
of  parameters  more  directly  related  to  the  quantities  obtained  in  direct  measure- 
ments; moreover,  we  find  that  this  new  set  of  independent  parameters  is  not 
explicitly  dependent  on  Te.  We  have  obtained  this  set  which  yields 


4-  = 4-  (5,  4>  , R,  t,,  S,,  9) 


(1.1^) 


This  expression  is  equivalent  to  1.13,  yet  it  exclusively  involves  dependences 
related  to  the  probe  and  to  ion  characteristics,  with  the  exception  of  6 which 
is  only  weakly  dependent  on  Te.  Using  this  appropriate  conversion  Sanmartin' s 
estimate  of  the  current  to  the  'finite'  cylindrical  probe,  when  subjected  to  the 
end  effect,  becomes 
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where 


T (e)  "-2—  a2  [X(p,  ,)] 
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(1.16) 


(1.17) 


a . ERF  ( J in  a2  J = 2/TV  R • (^  tt  + 0.6-rg2 


(1.18) 


We  have  elected  to  call  X( p. , tj ) the  peak-shape  function  since  it  exclusively 
describes  the  'nondimensional'  peak  structure.  The  expression  for  the  peak- 
shape  function  is  given  by 


with 


and 


X(p,tj)  h Xq(h,t,)  + | Xx(p,T)) 
X0(n,T|)  Xx(n,n)  axi 
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(1.20) 
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Since  no  general  analytical,  solution  exists  for  X(|i,tj),  it  has  to  he  evaluated 
numerically.  However  the  trivial  cases  |i  = 0 and  tj  = 0 yield  an  exact  analytical 
form  which  was  derived  by  Sanmartin  (Ref.  23) 


X(n  = OjTj)  = X(0,tj)  (1.22) 

X(n,tj  = 0)  = X(n50)  (1.23) 


with  the  following  limits 


Tj^O  -41  (cold  ions) 


(1.24) 


lim 

T)->  oo 


[X(0,Tj)] 
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(quasi- static) 

nj 


(1-25) 


1.1.3  Ion  Tengierature  Sensitivity 


The  marked  experimental  evidence 
with  the  typical  traces  shown  in  Fig.  1.3, 
narrow  peak  region  is  characterized  by  two 
height  of  the  peak  given  by 


of  the  end  effect  which  is  exemplified 
is  properly  described  by  1.13.  The 
important  peak  features,  namely  the 


height  = J(6  = 0)  s j^e  = 0)  • lf<y_  d2[x(n  = 0,T))]  (1.26) 

2R n/sT 
P 

and  full-width-at-half-maximum  (FWHM)  defined  such  that 

J ( 9 H \ J(e  = 0)  (I-27) 
or,  in  the  nondimen si onal  form 


X(n  = W.,tj)  = i X(p  = 0,tj) 
2 2 


(1.28) 
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such  that 


FWHM  = 


sm 


-l 


in  a2 


2(0^  - l) 


(1.29) 


The  peak- shape  function  X(p,rj)  describes  the  current  angular  response  by  the  way 
of  the  modified  angle  p,  and  the  current  sensitivity  to  the  ion  temperature  by 
way  of  tj.  The  overall  behaviour  of  X(p,Tj)  is  described  in  Fig.  1.4  for  the 
cold-ion  limit 


TJ  ->  0 or  S . — > °o 

1 1 

and  for  three  other  values  of  tj.  An  interesting  relation  can  be  found  by 
observing  that  all  the  peak-shape  functions  merge  with  the  cold-ion  case  for 
values  of  p smaller  than  their  half-width  values  pi , such  that 

2 

X(p\,  tj  f 0)  = X(pi,  cold) 

2 2 

Therefore,  using  Eq.  1.28,  we  have 

X(pi,  0)  = \ X(0,  tj)  (1.30) 

2 <L 

This  constitutes  a considerable  simplification  since  an  exact  analytical 
expression  for  I.l4  only  exists  for  the  two  trivial  cases  p = 0 and  tj  = 0.  The 
general  form  of  X(p,Tj)  must  in  fact  be  derived  numerically. 

Similarly  an  overall  description  of  the  peak- shape  function  can  be 
given  in  terms  of  the  various  angular  'landmarks’  which  characterize  the  current 
peak.  We  have  summarized  this  description  in  Fig.  1.5. 

The  sharp  current  response  obtained  by  the  finite  probe  about  the 
flow  axis  specifically  contains  the  ion  temperature  information  in  both  peak 
features.  Indeed,  once  the  physical  parameters  pertinent  to  the  probe  and  to 
the  ion  stream  are  found 


d,  i , u,  M,  n. , $ 

’ ’ ’ 1’  p 

we  have,  from  1.26  and  I. 29, 


J(6  = 0)  = J(5,  S.,  4>p) 


(1.31) 


01  = 0^(5,  S.  , 4>  ) 

22  1 P 


(1.32) 


The  sensitivity  of  these  peak  features  to  Si9  the  ion  speed  ratio, with  the 
geometric  aspect  ratio  as  a parameter,  can  be  developed  using  the  master  plots 
given  in  Figs.  1.7  and  1.8.  The  presentation  of  these  plots  is  simplified  by 
the  fact  that  the  nondimensional  peak  height  is,  to  first  approximation,  directly 
proportional  to  6 and  */$!  whereas  the  FWHM  is,  to  first  order,  independent  of 
5 and 

The  dependence  of  the  end  effect  on  can  be  obtained  in  terms  of  the 
parameter  Tg,  which  characterizes  the  true  finiteness  of  the  probe  in  terms  of 
the  flow  conditions.  The  nondimensional  current  to  the  parallel  cylindrical 
probe  obtained  from  1.26  is  plotted  in  Fig.  1.6  as  a function  of  t g for  the  range 
of  validity 


; 

f i 

I ' 

I > 


0 < r£  <3 


with  €*  as  a parameter.  This  parameter  describes  the  relative  size  of  the  probe 
radius  in  terms  of  a judicious  ratio  grouping  the  two  nondimensional  probe  lengths, 
namely  Tj  and  R.  We  get 


e*  = 


2R 


(1.33) 


(I.3M 

(1.35) 


This  last  quantity  is  a characteristic  length  which  is  defined  in  the  same 
fashion  as  the  Debye  length,  only  it  uses  the  bulk  ion  energy  instead  of  the 
electron  energy  since  the  electron  Debye  length  is  given  by 


e kT 

_o e 

e2  n. 

l 


Using  Fig.  1.7  one  verifies  that  j(0  = 0)  grows  with  S^,  whereas  the  current 
peak  decreases  with  increasing  e*  and  therefore  probe  radius,  under  similar  flow 
conditions.  Furthermore  the  end  effect  can  lead  to  a significant  departure  from 
the  'infinite'  probe  response  (j/j  (0)  > 10),  and  more  noticably  so  as  the  non- 
dimensional  length  of  the  probe  (Tg)  is  varied  from  small  values,  say  0.01  to 
values  near  unity.  For  values  of  Tg  beyond  1 the  effect  is  relatively  insensitive 
to  a change  of  the  true  probe  length,  that  is  to  changes  in  &,  u or  n^.  These 
observations  have  been  verified  in  Chapter  IV  where  variations  of  n^  were  accom- 
plished by  moving  axially  downstream  of  an  ion  source,  in  order  to  scan  values 
of  Tg  smaller  than  3.  Beyond  this  point  the  theory  reported  here  is  not 
appropriate  (see  Fig.  1.9(b))..  These  graphs  have  been  prepared  for  arbitrary 
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values  of  6 and.  C;.  which  are  typical  of  the  experimental  range  covered,  namely 
0.25  and  0.1  respectively,  and  for  a typical  flow  condition  for  the  ions  yielding 
the  scaling  length  u 1 with  a value  of  0.1  cm  1 . _Thus  the  current  to  the 
parallel  probe  can  be  scaled  directly  using  a factor  6 's^^/0.25*0.1,  whereas  no 
adjustment  is  necessary  on  Fig.  1.8,  provided  these  values  are  within  a range 
that  guarantees  the  validity  of  the  analytical  expressions.  Considerable  informa- 
tion is  also  extracted  from  these  plots  by  looking  at  the  individual  performance 
of  various  probes.  A series  of  probe  aspect  ratios  ($/ d)  are  presented  for 
comparison  for  a few  standard  diameters.  The  arrowheads  simply  locate  the  speed 
ratios  beyond  which  the  end  effect  will  convey  no  more  sensitivity  to  speed  ratio 
either  from  the  peak  height  or  the  peak  width.  Using  Fig.  1.7  one  can  also  obtain 
an  idea  of  the  performance  of  one  probe  geometry  relative  to  another.  In  partic- 
ular one  can  conclude  that  the  smallest  aspect  ratio  shown  (i.e.  22),  for  example, 
does  not  lead  to  a significant  end  effect  (j/j  (0)  ~ l)  and  consequently  would 
not  be  a good  choice  as  a diagnostic  tool.  From  these  plots  one  would  conclude 
that  a probe  with  the  largest  practical  aspect  ratio  would  perform  best.  However, 
since  the  interpretation  of  the  data  requires  the  value  of  the  unknown  quantity 
5 we  will  show  later  that  the  comparison  of  the  peak  height  and  FWHM  can  yield  a 
useful  approximation  of  this  value.  The  accurate  measurement  of  very  small  angular 
widths  then  sets  a practical  limitation.  In  terms  of  this  compromise  we  can 

conclude  here  that  aspect  ratio.:  in  the  neighbourhood  of  200  yield  quite  satisfac- 

tory results. 

1.1.4  Range  of  Validity 

In  order  to  ascribe  a limit  to  the  validity  of  the  scheme  reported  by 

Sanraartin  one  must  go  back  to  the  assumption  restricting  the  estimate  of  the 

electric  field  in  the  vicinity  of  the  probe  which  was  obtained  in  I. 10.  This 
time-dependent  result  was  actually  obtained  by  introducing  time  averages  for 
5(t)  and  /^(t)  for  a particular  region 

r 


. m _ 
p < — = p 
- r m 
P 


Therefore  if  we  call  tm  the  typical  drift-time,  along  the  probe,  of  ions  which 
■were  at  this  boundary  pm  at  t = 0,  then  the  estimates  remain  valid  for  the  time 
interval 


0 < t < t 

— — m 

which  Sanmartin  reports  to  be  equivalent  to 

0 < t < 3a"* 


i.e.  smaller  than  about  three  ion  plasma  periods.  Since  the  drift  time  t is 
simply  ( steady- flow) 


u 
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the  limits  to  the  range  of  validity  can  he  expressed  in  terns  of  the  nondimen si onal 
prohe  length 


0 1 rji  = (r  ) fipi  < 3 (1.36) 

Moreover  as_the  probe  is  turned  by  an  angle  0 with  respect  to  the  bulk  velocity 
direction  (u),  more  and  more  ions  which  would  have  been  collected  at  t < tm  are 
lost  because  of  their  increasing  transverse  momentum.  Sanmartin  demonstrates 
that  there  is  a limiting  angle  0pim  beyond  which  the  end  effect  will  disappear 
and  which  is  given  by  the  relation 


1 
2l2 


u sin(0iim)  - 2i>±2  - ci 2r 

where  Vj_  is  the  radial  velocity  of  the  ion  due  to  the  'modified'  potential  (60p) 
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The  limit  also  reads 
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Using  Sanmartin' s estimate  for  pm,  which  represents  the  region  of  the  sheath 
effectively  contributing  to  j 


„ 2R  b ^ 

P = — a b <5 

m 7 p 


(1.38) 


we  obtain  the  pertinent  range  of  validity  of  the  end  effect  theory  in  terms  of 
the  length  of  the  'finite'  probe  and  in  terms  of  the  angular  span  of  the  narrow 
current  peak: 


1 < £_  < 22  Jl 

- r - T„  7 


°<r£<  3 
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For  the  range  9^m  < 9 < 7r/2,  the  probe  collects  the  same  current  as  if  its  > 
length  were  part  of  an  infinitely  long  cylinder  or,  according  to  1.4 


j 


- d.(e> 


For  probe  lengths  t £ in  excess  of  3 the  probe  current  density  is  described 
accurately  by  Bettinger  and  Chen's  analysis,  provided  the  probe  is  parallel  to 
u.  It  is  therefore  important  to  note  the  considerable  refinement  which  is  brought 
by  Sanmartin's  analysis  to  the  description  of  the  end  effect  when  conqpared  with 
the  work  of  Bettinger  and  Chen.  It  is  reflected  in  the  appropriate  choice  of  an 
effective  region  of  the  sheath  which  contributes  to  j,  instead  of  the  whole  extent 
of  the  sheath.  This,  we  shall  see,  leads  to  the  correct  description  of  j for 
Tg  < 3.  However  B-C's  current  at  9 = 0 grows  indefinitely  (Fig.  1.9(b)).  Never- 
theless their  result  is  correct  for  the  asymptotic  limit  Tg  ->  <»,  since  by  defini- 
tion we  must  have  for  all  9 


lim 

rra 


^-(9) 

Joo 


We  also  notice  that  the  values  of  the  peak  half-widths  disagree  for  Tg  < 3 in 
Fig.  1.9(a)  where  we  show  that  1.32  approaches  the  value 


for  t.  ->  00 


Finally  this  refinement  has  led  to  the  description  of  a modified  field  about  the 
probe.  This  information  is  conveyed  by  5 which  is  the  adjusted  field-scale,  or 
characteristic  field  length,  and  which  was  approximated  by  Sanmartin: 


1 

(1.39) 


The  proper  evaluation  of  the  value  of  6 has  in  fact  proven  to  play  a vital  role 
in  the  foundations  of  the  method  to  measure  the  ion  speed  ratio  using  the  end 
effect.  This  method  will  be  elaborated  in  our  Chapter  IV. 


in  \fp~  - i 
m 2 


1.2  Some  Ion  Temperature  Measuring  Concepts 
1.2.1  Brief  Survey 

Since  this  study  is  reporting  a novel  method  to  measure  the  ion 
temperature  it  was  deemed  necessary  to  give  at  least  a schematic  view  of  some 
of  the  various  methods  which  were  already  available  to  probe  users.  However 
we  do  not  intend  to  consider  this  section  as  an  exhaustive  survey  of  the  diag- 
nostic methods  used  in  static  or  dynamic  plasmas.  In  fact,  a number  of  surveys 
have  been  carried  out  on  this  subject  (Refs.  34-4l),  in  particular,  the  most 
recent  one  conducted  by  Talbot  et  al  (Ref.  33). 

Most  methods  attempting  to  measure  T^  were  actually 'developed  in 
conjunction  with  possible  applications  to  the  sounding  of  the  closest  natural 
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plasma  source  at  hand  - namely,  the  earth’s  ionosphere.  Ground  based  schemes 
involving  incoherent  backs cattering  radar- sounding  and  spectroscopic  methods 
have  actually  shown  poor  resolution  and  irregular  reliability.  Though  in  situ 
measuring  schemes  have  been  quite  limited  geographically  they  constitute  (Refs. 
34-40)  the  bulk  of  the  experimental  efforts  attempted  on-board  sounding  rockets, 
balloons,  and/or  satellites.  Since  data  and  a fortiori  comparisons  between 
the  methods  are  much  less  abundant,  only  a few  candidates  are  presented  here. 

They  actually  fall  under  three  basic  categories:  cup  analyzers,  active  metallic 

probes  and  stimulated  interaction  methods. 

The  cup  or  ion- trap  analyzers  have  been  used  quite  'universally'  after 
the  theory  of  Whipple  (1959) • Though  they  constitute  simple  systems  they  lack 
the  accuracy  to  separate  multiple  ion  species  and  more  so,  their  use  is  somewhat 
doubtful  as  they  likely  introduce  serious  disturbances  into  the  plasma  ahead  of 
the  probe  (Ref.  4l) . A collimated  Faraday  cup  system  (Ref.  25)  has  recently  been 
proposed  as  an  excellent  diagnostic  tool  to  determine  the  ion  temperature  in  a 
flowing  plasma  (see  Appendix  D),  yet  it  represents  a more  refined  variation  of 
the  orifice  probes  (kinetic-theory-probes)  quite  often  flown  on  ionospheric 
sounding  packages  (Ref.  36). 

Some  more  recent  methods  not  developed  specifically  for  ionospheric 
probing  have  demonstrated,  at  least  analytically,  a promising  ability  to  obtain 
Tj_  under  various  plasma  regimes.  A hybrid  design  by  Laframboise  and  Parker 
(Ref.  42)  for  orbit-limited  current  collection  resembles  in  construction  multi- 
electrode geometries  studied  by  Medicus  (Ref.  43).  An  AC  (RF)  scheme  developed 
by  Rubinstein  and  Laframboise  (Ref.  44)  to  measure  Tj_  actually  presents  the 
noteworthy  advantage  of  being  less  prone  to  severe  mechanical  damage,  as  for  in- 
stance for  the  diagnostics  of  fusion  plasmas,  because  the  probe  is  made  repulsive 
for  each  charged  species  in  turn  and  is  collecting  only  small  'sample'  currents. 
Another  novel  scheme  involves  probe  measurements  in  strong  magnetically  confined 
plasmas  (Ref.  45).  Other  studies  have  been  undertaken  to  bridge  the  gap  widened 
by  non-collisionless  flow  problems:  two  of  these  Ti  measuring  candidates  in  more 

dense  plasmas  include  the  works  by  Shih  and  Levi  (Ref.  46)  and  that  of  Chubb  (MHD) 
(Ref.  47).  Finally,  the  scheme  reported  herein,  which  concerns  the  now  standard 
cylindrical  geometry  for  the  collisionless  regime,  where  an  anomalous  end  effect 
is  brought  to  fruition  in  determining  the  ion  temperature  in  a high  speed  plasma 
(Ref.  22). 


As  for  typical  proponents  of  the  last  category  we  can  refer  the  reader 
to  two  methods  used  extensively  in  this  laboratory  for  ionospheric  sounding  (Ref. 
48) . An  electron-beam  fluorescence  probe  (EBFP)  is  used  to  measure  the  rotational 
temperature  of  the  heavy  ionosphere  components  by  analyzing  the  fluorescent  back- 
scattered  light  stimulated  locally  by  an  electron  beam.  A second  method  is  an 
aerospectrometer  which  also  uses  local  stimulation  of  the  neutral  constituents. 

A nude  ionization- gauge  type  hot  filament- slit-biased  collector  system  is  made 
velocity  selective  (attitude  variable)  to  obtain  both  density  and  speed  ratio 
measurements  of  the  collected  species. 


1.2.2  The  Merits  of  the  End  Effect 


It  is  quite  noteworthy  in  fact  to  observe  that  direct  methods  to 
determine  the  ion  temperature  in  static  plasmas  are  inherently  scarce.  Further- 
more there  are  only  a few  of  them  which  are  relatively  successful  in  flow  cases. 


The  reason  for  this  scarceness  in  static  plasmas  is  easily  recognized  since  the 
region  of  the  probe  collection  characteristics  which  is  exclusively  sensitive 
to  Ti  is  , in  most  regimes,  almost  nonexistent.  The  significantly  less  mobile 
ions  only  produce  marginal  currents  relative  to  those  given  by  electrons  unless 
large  negative  biases  axe  used.  Yet  at  moderately  large  negative  biases  the 
stable  sheath  which  is  created  around  the  probe  forces  an  ion  to  lose  its 
autonomous  temperature.  In  fact  the  ion  must  go  through  the  pre- sheath  region 
thereby  accelerating,  up  to  the  hybrid  speed 

£F~ 

vi^  J-w 


according  to  Bohm's  criterion  for  a stable  sheath.  The  current  collection  then 
becomes  sensitive  to  electron  temperature. 

In  flow  cases  the  'odds'  are  actually  improved  because  the  ions  may 
already  have  a sufficiently  large  directed  energy  toward  the  sheath  to  satisfy 
Bohm's  criterion,  and  to  reach  the  probe.  Under  ideal  collisionless  conditions 
and  at  moderate  biases  they  do  so  with  basically  their  'virgin'  temperature. 
Nevertheless  the  sensitivity  of  the  measuring  device  must  allow  a refinement 
of  the  measurement  of  T^  with  a satisfactory  accuracy.  In  most  cases  this  last 
requirement  inevitably  calls  for  considerable  sophistication. 

The  physical  reason  why  the  end  effect  shows  a strong  ion  temperature 
sensitivity  can  be  summarized  as  follows.  The  ions  that  contribute  to  the  peak 
gain  an  easy  entrance,  in  the  thick  sheath  collisionless  regime,  through  the 
front  of  the  sheath  when  the  probe  is  parallel  to  the  stream.  Consequently, 
their  transverse  angular  momentum  is  largely  undisturbed  and  since  their 
collection  strongly  depends  on  the  transverse  momentum,  a marked  ion  temperature 
dependence  results. 

The  practical  importance  of  the  sensitivity  of  the  end  effect  to  T^  is 
clear.  It  offers  a simple  extension  of  the  use  of  cylindrical  Langmuir  probes. 

An  angular  scan  of  the  ion  collection  provides  simultaneously  the  ion  speed  ratio 
and  the  flow  direction. 


II.  DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 


II .1  The  Testing  Facility 

The  vacuum  testing  facility  used  for  the  entire  research  program 
reported  here  is  the  UTIAS  Low  Density  Plasma  Tunnel  (Refs.  53-5^)*  The  chamber 
can  be  operated  in  different  pumping  modes  from  the  high  mass-flows  of  the  order 
of  104  ni/sec  at  pressures  of  the  order  of  10  (iHg,  to  flows  of  the  order  of  10 
|ji/sec  near  blank-off  pressure  of  about  5 x 10~6  torr  (see  Fig.  II-l).  The  test 
section  is  about  2m  in  diameter  and  its  volume  is  ~ 15,000i.  The  chamber  was 
equipped  with  a new  diffusion  pump  to  give  a higher  pumping  speed  in  the  low 
pressure  regime  and  LN2  traps  (lmr  area)  were  used  in  conjunction  with  a 
pressurized  LN2  dewar  filling  system  to  minimize  oil  backstreaming  towards  the 
probing  area  and  the  accelerator  structure.  A water-cooled  pyramidal  trap  is 
mounted  on  the  centre  line  and  can  be  positioned  downstream  to  provide  a collect- 
ing area  for  the  ion  beam.  An  angle  of  30°  was  selected  to  minimize  backstreaming 
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of  ions  and  backsputtering  of  copper  atoms  towards  the  test  area  (Fig.  II-2). 

A small  chamber  is  coupled  to  the  main  chamber  to  enclose  the  entire  accelerator 
system  when  it  is  withdrawn  from  the  tunnel.  A gate  valve  isolates  the  test 
chamber.  An  extension  mechanism  is  used  to  retrieve  and  to  extend  the  whole 
accelerator  unit  into  the  test  section  so  the  main  vacuum  is  not  broken  during 
maintenance. 

A three  degree- of- freedom  traverse  (XYZ)  coupled  to  position  indicators 
allows  scanning  of  the  greater  part  of  the  test  section.  The  flexibility  of  the 
traversing  rig  is  greatly  inproved  by  the  use  of  another  actuator  system  which 
adds  pitch  (±  105°)  and  yaw  (±  15°)  motions.  The  linearity  of  the  XYZ  mechanism 
is  within  1$  and  the  monitor  has  a resolution  of  l/lOO  of  an  inch.  Angles  can 
be  set  to  within  a few  tenths  of  a degree.  A special  anti-backlash  sprocket 
chain-drive  mechanism  is  used  to  achieve  smooth  and  positive  motion  (Fig.  II-3)  . 
All  actuating  motors  are  of  low  voltage,  synchronous  type  to  prevent  any  arcing 
troubles,  as  the  entire  system  can  be  engulfed  in  a low  density  plasma. 

Mounting  and  alignment  of  the  accelerator  unit  is  simplified  by  using 
locating  devices.  A guiding  rail  system  within  the  smaller  chamber  is  ccopletely 
attached  to  the  door;  the  use  of  quick-disconnect  electrical  leads  and  a self- 
aligning mechanism  insures  trouble-free  and  accurate  repositioning.  The 
separately  pumped  chamber  provides  independent  testing  without  requiring  the 
use  of  the  tunnel  thus  allowing  shorter  turn-around  pumping  times. 

A testing  bay  (Fig.  II-l)  allows  monitoring  and  operating  of  all 
functions  without  significant  displacement.  Feedthrough  ports  and  sighting 
Pyrex  ports  also  add  to  its  flexibility.  Pressure  measurements  are  done  with  a 
Bayard-Alpert  type  nude  ionization  gauge  calibrated  for  N2.  A correcting  factor 
of  about  0.75  was  used  throughout  for  the  reported  work  with  argon  (Refs.  55,56). 


II. 2 The  Ion  Source 


II. 2.1  Background 


The  use  of  ion  accelerators  and  plasma  sources  in  vacuum  facilities 
of  this  type  has  been  applied  to  various  studies  ranging  from  ionospheric 
simulation  to  plasma-flow  interactions  of  many  kinds  (Refs.  19,57-59).  More 
particularly,  the  type  of  accelerators  developed  in  conjunction  with  ion  propul- 
sion technology  since  the  mid-sixties  has  been  a successful  candidate  in  the 
majority  of  these  cases.  The  ion  accelerator  used  for  this  study  was  developed 
in  this  laboratory  to  serve  as  a bench  prototype  for  an  ion  thruster  program 
initiated  in  1970-1971. 

The  accelerator  produces  a synthetic  plasma  of  high  velocity  when  a 
neutralizing  electron  current  is  injected  into  the  positive  ion  beam.  A schematic 
diagram  shown  in  Fig.  II-4  outlines  the  main  conqponents:  the  neutral  gas  source, 

the  plasma  source,  the  accelerating  and  beam-shaping  structure  and  the  neutralizer. 
The  present  design  benefited  from  the  various  scaling  relationships  available 
from  the  literature  (Ref.  60) . Two  nominal  beam  diameters  of  2.54  and  5.08  cm 
were  used  in  the  experiments.  The  smaller  diameter  was  obtained  by  masking  part 
of  the  source. 
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II.2.2  Gas  Source 


The  neutral  gas  source  is  controlled  by  a double  needle-valve  and  the 
mass-flow  is  monitored  by  noting  the  pressure  drop  with  time  in  a plenum  of  known 
volume.  Various  gases  have  been  used  successfully  (He,  N2>  Ar)  with  this  setup. 
Rare  gases  (Ar,  Xe,  He)  are  of  course  preferred  for  simplicity  of  operation. 

The  experiments  reported  here  were  restricted  to  argon  because  it  is  cheap  and  a 
wealth  of  background  data  is  available.  The  stability  of  the  neutral  flow  is 
excellent.  The  method  can  be  very  accurate  and  provides  very  low  mass-flow  rates 
down  to  1.5  ni/sec  for  argon.  A distributor  plenum  is  attached  to  the  ion  source 
just  upstream  of  the  plasma  chamber  to  increase  the  homogeneity  of  the  flow  (see 
Appendix  A) . 


II. 2. 3 Discharge  Chamber 

The  gas  enters  the  plasma  discharge  chamber  evenly  and  is  ionized  by 
electron  bombardment  in  a hot-cathode  reflex  discharge  (Fig.  II-7)  . A modest 
magnetic  field  (~  50  Gauss)  is  used  to  improve  the  ionization  efficiency  of  the 
hot  cathode- electrons  in  a crossed- field  configuration  such  that  their  total 
effective  distance  travelled  is  increased  manyfold  and  the  radial  drift 
component  minimized.  The  anode-cathode  voltage  drcp  is  optimized  to  satisfy  the 
main  requirement  that  the  discharge  voltage,  Vp  be  at  least  twice  the  value  of 
the  ionization  potential  (~  15-7  eV  for  first  ionization  of  argon)  and  is 
nominally  set  at  50  volts  which  corresponds  to  the  optimum  energy  for  argon 
(Ref.  79) • The  directly  heated  filament  is  positioned  within  a small  cylindrical 
pole-piece  which  enhances  the  magnetic  field  near  the  filament  to  minimize  elec- 
tron losses  to  the  anode.  Many  filament  configurations  were  considered,  namely, 
hairpin,  coil,  button,  but  the  coiled  type  was  preferred.  A quick-access  filament 
support  with  two  filaments  mounted  at  all  times  was  used.  Typical  lifetimes  in 
argon  are  about  30  hours  (Fig.  II-8)  at  emission  densities  of  about  0.35  anq?/cm2 
(0. 008-0 .010  in.  tungsten  and  0.007  in.  tantalum).  The  entire  internal  surface 
of  the  chamber  including  the  perforated  end-plates  are  mirror  polished  to  minimize 
heat  losses.  As  the  electron  emission  (Richardson-Dushman  law)  grows  faster  with 
temperature  than  pure  radiative  power  (ideal  black-body),  satisfactory  operation 
of  discharges  of  this  type  normally  do  not  require  cooling  of  the  discharge 
chamber  walls,  provided  of  course  that  the  materials  are  carefully  selected  for 
all  components.  Temperatures  slightly  in  excess  of  700°K  have  been  measured  on 
the  discharge  anode  structure. 


II  .2.4  Accelerator 


The  accelerating  unit  consists  of  a flat-plate  electrostatic  system. 

Two  perforated  match-drilled  plates  contain  a regular  pattern  of  holes  in  hexa- 
gonal symmetry  (Fig.  II-8)  . The  required  precision  necessitated  a jig  production 
scheme.  The  selection  of  the  hole  sizes  grid  spacing  and  plate  thickness  is 
critical,  to  the  proper  operation  of  the  accelerator  and  will  be  dealt  with  under 
the  flow-field  considerations.  The  nominal  sizes  were:  1.25  mm  thickness 

(stainless  steel),  1.56  cm  diameter  holes  and  a 1.0  mm  spacing.  A mask  with  a 
2.54  cm  diameter  opening  can  be  inserted  on  the  discharge  side  of  the  first  grid 
to  provide  the  smaller  beam  diameter  (Fig.  II-8) . A proper  designing  of  all  the 
accelerator  parts  is  also  required  for  high-voltage  breakdown  (Paschen)  protection. 
A small  cylindrical  grounded  shield  about  twice  the  diameter  of  the  accelerator, 
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protruding  into  the  main  vacuum  chamber,  was  found  successful  in  providing 
reduction  of  Paschen  breakdowns.  The  electrical  insulating  materials  used 
include  Teflon,  lava,  Ajiz03,  mica,  mylar,  Kapton,  Pyrex  and  porcelain  (ceramic) 
cements. 

•» 

II.  2. 5 Neutralizer 

The  neutralizing  scheme  adopted  was  a directly  heated  filament 
immersed  at  the  beam  edge.  Straight  wires,  strips  and  coiled  wires  have  been 
tested  and  the  latter  were  preferred.  Early  studies  used  a hot-filament  array 
within  the  beam  downstream  of  the  source  in  conjunction  with  a screen-cage 
enclosing  the  whole  beam.  The  method  was  judged  too  awkward  since  it  limited 
access  to  the  full  test- section  and  precluded  the  probing  of  the  entire  beam 
envelope.  The  preferred  neutralizer  filament  was  mounted  directly  on  the 
accelerator  structure  as  it  allowed  maintenance  without  breaking  the  main 
chamber  vacuum  and  provided  trouble-free  operation.  The  measured  neutralizing 
current  was  usually  in  excess  of  85 % of  but  typically  less  than  the  emitted  ion 
current.  It  is  thought,  as  reported  by  other  investigators  that  a small  per- 
centage of  background  electrons  is  injected  into  the  positive  beam  from  nearby 
ground  surfaces,  from  secondary  electron  sputtering  and  other  collision  mechanisms 
within  the  beam. 


11.2.6  Operating  Conditions 

The  selection  of  the  basic  operating  parameters  (Figs.  II-6,  II-9)  was 
made  with  prime  concern  to  obtain  a reliable  and  reproducible  range  of  experimental 
conditions.  The  range  of  beam  current  available  goes,  with  the  nominal  plate 
spacing,  roughly  from  0.6  to  12.5  mA;  a standard  setting  of  2 mA  was  taken.  It 
allowed  stable  operation,  in  conjunction  with  a filament  current  regulator,  at  a 
point  minimizing  neutral  gas  losses,  allowing  extended  running  time  with  minimum 
attention  and  it  also  provided  a reference  current  compatible  with  a wide  range 
of  probe  currents  for  beam  studies.  Care  was  also  taken  to  operate  within  'safe' 
boundaries  with  regard  to  excessive  noise  on  the  ion  stream.  The  nominal  dis- 
charge was  operated  at  50  volts  (50-100  mA.)  with  an  electromagnet  current  of 
5-6  anps.  The  mass-flow  was  typically  set  around  9*5  equivalent  milli-anps  of 
U.H.P.  argon*. 

As  we  shall  see  shortly  the  beam  homogeneity,  its  stability,  its 
symmetry  along  with  its  geometrical  features  in  general,  are  of  paramount 
importance  to  the  study  described  here.  The  normal  high-voltage  configuration 
(Fig.  II-4)  was  used  extensively  for  this  experiment  (150-1200  eV  ions).  A 
low-voltage  configuration  (10-300  eV)  using  a single-plate  extractor  was  also 
used  (Fig.  II-5) . Some  qualitative  impressions  can  be  obtained  by  actual  beam 
pictures.  Photographic  records  of  beam  envelopes  are  presented  for  comparison 
between  the  different  modes  of  operation  (Fig.  11-10) . Severe  arcing  is  not 
common . However  intermittent  and  unpredictable  occurrence  is  usually  linked  to 
extended  operation  without  maintenance.  A shadow- shield  was  added  to  protect 
the  insulators  from  neutralizer  material  vaporization. 

An  interesting  feature  of  severely  un-neutralized  beams  or  arc-breakdown 
is  the  fanning  out  of  the  beam  (l80°)  often  accompanied  by  arcing  at  the  tunnel 

* 1 equivalent  mA  (argon)  ~ 0.l8  pil/sec  or  1.5  mgm/hour. 
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back-plate  as  for  the  case  shown  in  Fig.  II-10(a);  such  occurrences  only  take 
place  when  the  neutralizer  is  turned  off  completely,  as  is  the  case  when 
taking  photographic  records,  for  the  filament  provides  a natural  30  to  50 
watt  incandescent  'bulb' . 


II. 3 Langmuir  Probes 


i 

( 


II. 3.1  Construction 

The  use  of  Langmuir  probes  in  plasma  studies  is  well  documented 
(see  Chapter  I and  Refs.  l-l6).  However,  there  does  not  seem  to  be  a single 
best  way  for  the  construction  and  for  the  operational  'rules  of  the  game' 
since  so  many  various  applications  and  requirements  exist.  Many  practical 
experimental  difficulties  such  as  probe  contamination  are  still  debated.  We 
refer  to  Appendix  B which  explicitly  considers  a few  anomalous  effects  of 
some  concern  to  us  here.  One  very  practical  difficulty  involves  the  construc- 
tion of  high  aspect  ratio  {&/ d)  probes  such  as  the  ones  shown  in  Fig.  II-3. 

We  have  presented  in  Fig.  Il-ll(a)  a list  of  the  various  probes  used  in 
conjunction  with  this  experimental  investigation,  and  in  Fig.  Il-ll(b)  an 
outline  of  the  many  experimental  runs  recorded  with  these  probes,  in  terms 
of  the  different  source  parameters  and  probe  parameters  which  were  required. 
Figure  11-12  shows  the  basic  design  features  of  the  probes  used  for  this 
investigation  which  have  been  standardized  for  e ase  of  construction  and 
mounting.  These  probes  are  rugged  and  can  sustain  high  temperatures  yet 
their  handling  and  mounting  requires  caution,  especially  with  the  most  slender 
ones.  The  mounting  on  the  scanning  rake  (Fig.  II-3)  requires  the  positioning 
of  the  probe  to  within  approximately  half  a diameter  with  respect  to  the 
centre-line.  The  probe  dimensions  are  measured  to  within  2$  of  their  length 
and  5$  of  their  diameter. 

Satisfactory  straightness  of  the  thin  probes  is  assured  by  using  a 
relatively  effective  straightening  method  which  was  suggested  in  Ref.  6l.  A 
dead  weight  is  suspended  from  a length  of  tungsten  or  tantalum  wire  and 
resistive  or  flame-heating  is  used  for  a short  length  of  time  to  provide  the 
straightening  effect.  A scale  is  removed  afterward  using  the  electrolytic 
action  of  a NaOH  cleaning  solution  with  a 1 to  2 volts  voltage  difference 
for  a few  seconds.  The  surface  is  then  degreased  and  handling  is  restricted. 
Flexures  or  'kinks'  of  less  than  one  diameter  can  be  achieved  with  the  smaller 
size  wires  (0.01  cm).  However,  since  the  natural  bending  of  the  'cantilevered' 
cylindrical  structure  which  the  conducting  wire  constitutes  can  be  in  excess 
of  one  diameter  for  the  long  slim  cases,  this  has  been  considered.  Indeed, 
as  we  demonstrate  in  Appendix  H,  since  the  nondimensional  deflection  is  propor- 
tional to  the  aspect-ratio  (R  = £/d)  to  the  third-power,  a deflection  in  excess 
of  5 diameters  is  expected  for  an  ideally  'straight'  horizontal  wire  at  R = 472 
for  d ~ 0.01  cm  (0 .004  in)  . This  leads  to  an  angular  error  at  the  tip  of 
nearly  0.9° . 


II. 3 .2  Probe  Cleaning 


Further  probe  cleaning  is  required  after  the  fairly  polished  finish 
is  achieved  using  base  baths  (NaOH),  in  order  to  ensure  repeatable  and 
believable  probe  measurements.  Pure  heating  is  insufficient  and  usually 
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impracticable  at  the  densities  involved  (n^  < 10®  cm  3) . Resistive  heating  to 
about  250°K  for  10  minutes  (~  10  5 Kw/cm2)  has  been  reported  in  Ref.  62  by 
Pigache  to  be  effective  at  densities  of  the  order  of  10 6 cm  3 . However,  as  is 
reported  by  the  same  author,  ion-bombardment  cleaning  is  more  practical  and 
furthermore  readily  available  by  using  the  ion  beam  from  the  accelerator  itself. 

Pigache  reports  2 KeV  argon-ion  cleaning  for  the  same  period  of  time  to  be 
equally  effective.  Though  his  results  were  for  stainless  steel,  we  estimate 
that  for  tungsten  at  the  densities  involved  in  this  present  experiment  a 
cleaning  time  of  ten  minutes  would  require  ^ 1 KeV  argon  current  density  of 
80  pA/cm2  giving  a heating  power  of  8 x 1C>  kW/cm2.  Usually  this  power  level 
was  exceeded  by  biasing  the  probe  by  about  450  volts  negative  at  a beam  energy 
of  1 KeV  such  that  5 to  6 minutes  were  deemed  sufficient.  Indeed,  a procedure 
of  regular  cleaning  and  keeping  the  probe  continuously  within  the  beam  gave 
satisfactory  results.  No  random  current  jumping  or  excessive  slow  drifting  of 
the  measurements  was  noted  for  periods  in  excess  of  20  minutes  after  ion- 
bombardment  cleaning.  A simple  model  for  probe  cleaning  was  attempted  in 
Appendix  C. 

II. 4 Cup  Probes 

II. 4.1  Button  Probes  I 

Simple  flat-button  probes  at  the  bottom  of  insulating  cups  (see  Figs.  k 

II-3  and  II-l4(b))  were  constructed  to  measure  the  slow  background  ions  present 
within  the  beam.  These  probes  can  be  electrically  biased.  They  are  mounted  on 
the  same  traversing- scanning  rig  as  the  cylindrical  wire  probes. 


11.4.2  Collimated  Cups 

The  schematic  of  a collimated  Faraday  cup  given  in  Fig.  A-2  shows  the 
orifice  geometry  selected  to  sight  the  smallest  area  possible  at  the  accelerator 
exit-plane  for  beamlet*resolution,  yet  large  enough  to  yield  a reasonable  current 
level,  typically  of  1 n amp.  To  keep  the  dimensions  small  a compromise  was 
arrived  at  of  a 1.75°  full-angle  aperture.  The  orifice  size  is  0.08  cm  and  the 
second  orifice  can  be  biased  to  provide  a repelling  potential  for  the  secondary 
electrons  originating  from  the  collimator.  The  collector  was  normally  set 
floating  and  the  'grid'  biased  at  -9  volts.  An  adapter  is  used  to  mount  the 
miniature  cup  on  the  standard  rake.  The  cups  have  been  used  to  define  the  beam 
edge  and  to  study  the  spreading  of  the  beam.  It  also  provides  an  interesting 
diagnostic  method  to  determine  the  local  flow  velocity  direction  and  possibly 
the  ion  speed  ratio  (see  Chapter  III,  Section  III. 1.3  and  Appendix  D) . 
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II. 5 Instrumentation 


II. 5.1  Ancillary  Equipment 

The  operation  of  the  pumping  facility  is  accomplished  by  a fully 
interlocked  fail-safe  control  system.  The  XYZ  traversing  mechanism  can  be 
operated  using  a hand-held  multi -function  control  box.  The  scanning  mechanism 
(Fig.  II-13(b))  provides  both  digital  and  analogue  signals  to  be  monitored  at 
an  X-Y  recorder. 

*see  Figs.  III-l  and  III-2 . 

19 


11.5*2  Source  and  Accelerator 

A table  shown  in  Fig.  IIt  9 gives  the  basic  requirements  for  the  operation 
of  the  accelerator.  The  basic  feature  for  the  extended  stable  operation  of  the 
plasma  source  is  the  beam  current  regulating  system  based  on  a very  reliable 
design  used  by  the  UTIAS  Rocket  Group  for  their  flight-proved  electron  guns,  in 
which  the  discharge  filament  current  is  controlled  to  keep  the  beam  current 
constant.  The  circuit  is  protected  from  high-voltage  (Fig.  II-13(a))  surges 
and  is  current  limited.  The  filament  current  is  provided  by  a floating  12  volt 
automotive  battery  and  beam  current  regulation  can  be  achieved  over  a voltage 
range  of  about  120  volts  to  1800  volts  at  2 mA  ion  emission  for  extended  periods 
(>  3 hours).  A schematic  diagram  (Fig.  11-15)  describes  the  initial (warm- up) 
operation  of  the  source  before  the  beam  is  turned  on.  It  simply  involves  running 
the  regulator  with  an  external  'dummy'  calibrated  current  to  allow  the  discharge 
to  stabilize  before  the  high-voltage  is  turned  on. 


II. 5. 3 High-Voltage  Protection 

Safe  operation  with  positive  high-voltage  on  the  discharge  supplies  is 
allowed  by  the  use  of  plexiglass  to  insulate  the  supplies  and  their  controls  and 
by  the  use  of  a 1-to-l  10  KV  line  isolating  transformer.  A commercial  regulated 
power  supply  used  for  the  discharge  required  extra  care  since  an  arc-over  to 
ground  from  the  positive  high-voltage  supply  would  put  a transient  voltage  spike 
on  its  output  stage  and  destroy  it.  The  integrity  of  this  supply  was  preserved 
during  the  whole  extent  of  this  study  by  the  use  of  a simple  and  inexpensive 
surge-limiter  circuit  (Fig.  II-16) . A bank  of  gas-filled  surge  limiters  (V  ~ 70- 
90  volts)  is  also  put  across  the  output  terminals  to  limit  any  surges  beyond  the 
typical  operating  level  (50  volt)  with  a response  of  better  than  5 KV/psec. 

II. 5.*+  Probe  Circuitry 

The  single-Langmuir  probe  circuit  is  a staidard  balanced  bridge  design 
(Fig.  11-17)  with  a slow  manual  voltage  sweep  (X-axis).  The  measuring  leg  (Y- 
axis)  is  recording  better  than  99 lo  of  the  probe  current  signal  at  10  KG.  This 
circuit,  however,  contains  an  important  refinement  in  practice  to  provide  the 
ability  to  re- trim  the  bridge  frequently  when  low  currents  are  obtained  (typi- 
cally Z 10  6 amp).  Indeed,  as  the  input  impedance  of  the  Y preamp  on  the  chart 
recorder  can  vary  slightly  with  scale  factor  and  since  large  bias  voltages  are 
utilized,  small  imbalances  can  then  be  readily  corrected  to  avoid  significant 
offset  errors.  The  frequency  response  of  this  entire  circuit  is  limited  at 
1 KHz.  A 5 hF  capacitor  across  the  Y-axis  terminals  of  the  recorder  was 
occasionally  necessary  to  reduce  noise  at  low  levels. 

For  measurements  at  constant  voltage  a sensitive  current- to-voltage 
converter  was  preferred.  It  allowed  the  measurement  of  currents  of  10"9  to  10-3 
Anp  at  1 volt  full-scale  output  with  negligible  drift  for  extended  periods.  The 
circuit  response  is  kept  below  100  Hz  (typically  30  Hz)  to  minimize  high-frequency 
noise  when  used  in  conjunction  with  the  probe  pitch  and  yaw  scans  (Fig.  II-l4(a)). 
This  frequency  response  is  quite  sufficient  to  monitor  the  signals  since  the  ± 105° 
scanning  is  accomplished  in  a relatively  long  period  - typically  a few  seconds. 
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Standard  current  and  voltage  monitoring  meters  are  used  with  accuracies 
of  better  than  3$  full-scale,  throughout.  The  XY  recorder  is  a series  7001A 
Hewlett-Packard  (metric)  unit  and  the  re- calibration  procedures  were  performed 
regularly. 

III.  FLOW-FIELD  CHARACTERISTICS 

111 .1  Ion  Temperature  Simulation 

111. 1.1  Ion  Velocity  Distribution 

The  correct  simulation  of  an  ideal  thermal-ion  stream  with  drifting 
velocity  (bulk  velocity)  u relative  to  some  object  to  be  tested  requires  that  the 
flow  be  described,  in  the  drifting  frame  of  reference,  by  a standard  Maxwellian 
distribution  function  giving  the  phase-space  density 
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In  the  laboratory  reference  system  (primed)  therefore, we  can  express  the  resulting 
velocity 
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such  that  the  'longitudinal'  thermalized  population  can  be  described  by 
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or  using  the  energy  relation  E = l/2  Mu2  and  E , = l/2  Mv 
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This  longitudinal  distribution  is  characterized  by  a thermal  energy  spread  which  is 
roughly  described  by  the  points  at  which 


(e;)2  - (Eq)2  = (kT.) 


(III. 5) 


If  we  assume  that  E0  » kTi?  i.e.  since  Sj_  » 1,  then  we  obtain  an  approximation 
for  this  full  width  between  the  (l/e)  points  (e  = exp[0.1]) 


or 
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(ill. 7) 


where  AW  is  the  natural  width  of  the  distribution  in  the  fixed  laboratory  frame 

(2^). 

The  perfect  simulation  of  the  ion  temperature  for  the  transverse  component 
however  is  simplified  since  there  is  no  dependence  on  the  translation  motion.  Thus 
it  implies  that 
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and  a full  width  identical  to  the  natural  energy  spread 


AW'  = AW  s 2KT. 
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(III. 10) 


Therefore  the  longitudinal  simulation  requires  a considerably  broader  energy  width 
than  the  transverse  case  since  in  most  applications  E0  » kT^  in  Eq.  (III. 7).  For 
instance,  to  simulate  the  correct  motion  of  a vehicle  at  satellite  velocity  (~  10 
km/sec)  through  a nitrogen  plasma  at  about  1200°K  (0.1  eV)  the  N2+  stream  in  the 
test  tunnel  should  show  an  energy  width  of  nearly  5 eV  centred  about  a main  energy 
of  15  eV  for  the  longitudinal  component,  and  a width. of  only  0.2  eV  for  the 
transverse  thermal  component. 


Ill .1.2  Plasma  Stream  Model 

The  relative  motion  between  a probe  and  the  ions  and  electrons  of  a 
streaming  plasma  is  the  result  of  the  directed,  or  bulk  velocity,  and  the  thermal 
velocity  of  the  particles  as  we  have  seen  above,  tore  specifically  the  plasmas 
of  interest  to  this  report  satisfy  the  mesothermal  condition 
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(III. 11) 


which  applies  to  a great  number  of  natural  and  laboratory  plasmas  and  in  particular 
to  the  ionosphere.  It  indicates  that  the  relative  velocity  between  the  probe  and 
the  ions  is  essentially  u with  a small  random  thermal  component  <v^>  added  to  it 
whereas  the  probe-electron  relative  velocity  is  totally  randomized  in  direction 
and  has  a magnitude  given  by  <Vg>.  This  also  means  that,  in  ord6r  to  achieve  an 
accurate  simulation,  the  small  random  velocities  represent  a refinement  which  may 
be  omitted  in  many  cases,  yet  it  cannot  be  disregarded  in  the  study  of  the  beam 
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structure  or,  as  is  the  case  here,  in  the  description  of  a probe  effect  which  is 
very  sensitive  to  exactly  this  transverse  random  velocity  component. 

However,  the  particular  difficulties  associated  with  the  generation  of 
an  extended  plasma  beam  with  ideal  Maxwellian  features,  at  drifting  velocities 
approaching  satellite  values,  are  such  that  a situation  short  of  the  ideal  has  to 
be  accepted  for  experimental  purposes . 

The  beam  produced  by  an  electron-bombardment- type  ion  thruster  (Chapter 
II)  actually  meets  most  of  the  requirements  of  simplicity  in  that  respect,  and  the 
distribution  function  can  be  modelled  satisfactorily  in  most  aspects.  The  electrons 
of  the  well  defined  plasma  stream  are  essentially  totally  randomized  at  a temperature  . 
close  to  the  neutralizing  wire  temperature  (see  Section  III .2.2) 

< v > ~ Ttt  ~ i eV 
e W 3 

A significant  part  of  the  longitudinal  ion  energy  spread  comes  from  the  fact  that 
the  ions  are  extracted  non-ideally  through  a complex  sheath  geometry  in  a region 
of  the  accelerator  where  large  potential  gradients texist . Furthermore,  this  spread 
is  affected  by  a complex  interaction  or  cross- coupling  between  the  mags- flow  and 
the  source  discharge  characteristics.  It  can  be  estimated  (Ref.  28)  that  this 
longitudinal  spread  is  typically  less  than  about  5%  of  the  main  accelerator  voltage 
or  about  ± 10  eV  for  an  accelerator  potential  of  500  volts.  According  to  Eq.  III. 7 
this  would  be  just  about  sufficient  to  properly  simulate  an  ion  temperature  of 
0.1  eV  at  a beam  energy  of  500  eV  (e.g.  50  km/sec  argon).  However  no  attempt  was 
made  to  purposely  tailor  this  longitudinal  spread,  as  it  is  expected  to  be  of  little 
consequence  in  this  investigation,  where  only  the  transverse  spread  is  of  paramount 
importance. 

An  exact  description  of  the  physical  mechanisms  involved  in  the  generation 
of  the  natural  transverse  energy  spread  of  this  same  beam  is  considerably  complex. 

The  main  beam  can  be  depicted  by  a superposition  of  a series  of  small 
discrete  beamlets  extracted  by  each  individual  circular  aperture  as  shown  in  Fig. 
III-2.  The  small  beamlets  are  expanding  with  a limiting  divergence  angle  2 y and 
are  merging  smoothly  into  a well  defined  beam  which  is  contained  within  an  angle 
2y. 


The  ion  contribution  to  a point  on  axis  at  large  distances  downstream 
is  defined  in  Fig.  III-l  for  a circular  aperture  As^  located  at  the  maximum  source 
radius  r0.  The  transverse  conqponent  of  the  velocity,  for  an  ion  released  from  As^ 
and  collected  at  an  axial  distance  z is  given  by 


•vkj_  = \ sine 


(III. 12) 


provided  one  assumes  collisionless  and  field-free  conditions  within  the  beam.  The 
very  nature  of  a cylindrically  symmetric  source  of  finite  extent  actually  provides 
a maximum  transverse  velocity  component  for  an  ion  arriving  at  a position  (z,o) 
under  these  conditions,  namely 
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where  u is  some  average  directed  velocity  for  the  ions. 


It  is  also  reasonable  to  assume  that  every  one  of  the  closely  packed 
apertures  contributes  to  the  main  beam  as  a separate  source  and  that  since  only 
small  angles  are  involved,  one  can  expect  that  an  ion  collected  at  a large  axial 
distance  (in  terms  of  the  source  radius  r0)  can  originate  from  the  separate  source 
with  an  equal  probability.  We  can  also  assume  that  the  velocity  distribution  from 
each  beamlet  will  be  zero  beyond  the  maximum  value  given  by  III. 13.  Therefore, 
following  Hester  and  Sonin  (Ref.  l8)  we  can  represent  the  distribution  for  the 
magnitude  for  the  transverse  velocity  in  the  x-direction  (by  analogy  the  same 
treatment  holds  for  the  y-direction)  by 


f(cx)dcx 


(III.14) 


To  relate  this  to  an  equivalent  Maxwellian,  we  make  the  observation  that 
the  average  magnitude  of  a velocity  component  equals  half  the  mean  random  speed, 

i .e. 


c 


therefore  calling  this  average  transverse  velocity  <Vj>  we  have 

c = <v,>  = 2~\c  I 
-L  1 x 1 


and 


c 

max 


v± 

max 


where 


v± 

max 


u 


u 


r 

o 
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We  get  the  average  transverse  velocity  of  the  ions  arriving  at  z 


< V|  > ~ •=-  


Jtt 


(III.16) 


An  effective  ion  temperature  is  then  defined  assuming  that  a synthetic  thermaliza- 
tion  of  the  ion  papulation  is  satisfactorily  achieved.  The  Maxwellian  average  is 
then  equaled  to  the  average  found 
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which  yields 

T ~ . f- °u")  (III. 17) 

X±  9k^  \ z / 

Similarly  a synthetic  ion  speed  ratio  is  defined  using  the  geometrical  model  with 
the  expression 


si  = ~ = ¥ ( T 
geo  c^  v o 


(III.18) 


For  example,  a 1 KeV  singly-ionized  argon  beam  would  simulate  a speed  ratio  of 
about  80  at  an  axial  distance  z/r0  ~ 40.  An  ion  temperature  of  about  0.5  eV 
would  be  simulated  in  the  transverse  plane  according  to  III .17. 


III. 1.3  Verification  of  the  Source-Flow  Model 

In  order  to  summarize  the  assumptions  underlying  the  beam  model,  they  can 
be  listed  briefly  as  follows: 

1.  The  flow  conditions  satisfy  the  collisionless  criteria  (see  Fig.  1-1  and 
Appendix  E) . 

2.  The  charged  species  are  streaming  in  field- free  space  away  from  the  extractor, 
assuming  that  no  residual  fields  are  present  and  that  neutralization  is 

complete. 

' 

3.  The  beam  density  is  approximated  using  a conical  expansion  model  with  a 
divergence  angle  y and  falls  off  roughly  as  the  inverse  of  the  axial  distance 

U/ror2. 

4.  The  geometrical  centre-line  is  preserved  throughout  (no  misalignment  and/or 
vectoring)  such  that  the  current  density  and  the  divergence  at  any  given 
radius  from  the  centre-line  are  axially  symmetrical. 

5.  The  angular  density  at  a given  radius  from  the  centre-line  is  symmetrical 
around  a contributing  solid  angle  sin0d8. 

6.  The  area  contributing  to  the  probe  current  is  large  enough  to  average  out 
single-hole  effects . 

7.  The  test  area  at  the  probe  is  small  enough  that  changes  in  current  density 
across  it  average  out  to  a value  representative  of  the  value  at  the  centre. 

The  conditions  to  validate  the  first  assumption  is  chosen  by  design  since  the 
collisionless  flow  conditions  are  pertinent  to  this  study.  These  conditions 
satisfy  the  criteria  discussed  in  our  Appendix  E.  The  restrictions  listed  under 
numbers  4 to  7 follow  from  the  symmetry  and  from  the  physical  dimensions  of  the 
source  relative  to  the  probe.  They  are  for  all  intents  and  purposes  straight- 
forward. Assumptions  2 and  3,  however,  imply  a considerable  combination  of 
physical  mechanisms;  they  are  the  ones  we  intend  to  discuss  more  specifically  in 
this  section. 


I 

t 
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An  ideal  description  of  the  beam  structure  is  depicted  in  Fig.  III-2. 

The  formation  of  the  beam  involves  the  interaction  of  the  large  field  gradients 
responsible  for  the  electrostatic  lens  action  of  the  apertures,  along  with  the 
space- charge  spreading  of  the  ion  beamlets  which  are  shaped  and  extracted  by  these 
fields  in  the  immediate  vicinity  of  the  grids.  Further  away,  the  ion  beam  is 
neutralized  by  a colony  of  thermal  electrons  emitted  by  an  incandescent  wire  in  a 
region  where  residual  fields  within  the  beam  are  small. 

The  divergence  angle  7 can  be  approximated  by  the  simplified  model  of 
Davisson-Calwick  (Ref.  77)  for  a pair  of  circular  grids  of  diameter  B,  and  where 
h and  s are  the  sheath  thickness  and  the  acceleration  distance,  respectively 


(III .19) 


The  angle  0^  is  the  angle  at  which  the  ions  are  assumed  to  be  released  from  the 
sheath,  with  basically  zero  energy.  The  parameter  a represents  an  important  ratio 
which  sets  the  optical  property  of  the  lens  in  terms  of  the  value  of  the  voltages 
V+  and  V on  each  grid 


a = 


V+  + |V_ 


Typically  we  have 


a ~ 0.7 


7 ~ 10°  - l4° 


which  implies  a value  of  h/s  ~ 0.04  if  we  assume  a plane  sheath  (9>!  = 0°).  Using 
some  nominal  source  estimates  of  the  Debye  length  from  values  for  che  electron 
temperature  and  density  within  the  source  (Te  ~ 3 eV  and  np  ~ 10 10  cm'3)  we  obtain 
a value  of  h/s  ~ 0.03  for  s ~ 4 mm.  This  result  suggests  that  the  divergence  model 
is  reasonably  accurate. 

In  general,  the  real  beam  shape  only  approximates  the  beam  geometry 
described  in  Fig.  III-2.  For  instance,  Fig.  II-10(e),(f)  shows  that  the  beam  edge 
is  not  truly  conical.  In  fact  the  beam  waist  is  slightly  constricted  near  the 
grids  and  gradually  expands.  This  is  resulting  from  the  convergent-divergent 
action  of  the  two  metal  grids  and  of  a virtual  electrode  situated  near  the  neu- 
tralizing point.  Only  a marginal  space-charge  effect  is  believed  to  be  present 
beyond  that  point  because  of  imperfect  neutralization.  This  can  be  demonstrated 
using  a fairly  accurate  model  for  the  space- charge  spreading  of  an  ideal  parallel 
ion  beam.  The  variable  of  interest  is  a product  of  the  beam  aspect  ratio  R , 
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D0/s  th^  ratio  of  source  diameter  and  acceleration  distance,  and  the  neutralization 
fraction  T (Ref.  69) 

01  Vl  - T (III. 20) 

Typically,  since  d{,  ~ 20  and  the  measured  divergence  angle  7 is  roughly  10*-  14°, 
as  indicated  in  Fig.  III-3,  the  corresponding  neutralization  would  be  about  99 *8$ 
complete  if  it  alone  were  to  account  for  this  divergence.  This  is  of  course  beyond 
our  measuring  accuracy;  nevertheless  it  makes  it  plausible  that  the  actual  curvature 
of  the  beam,  as  it  leaves  the  accelerator,  is  caused  by  these  straightforward  mech- 
anisms . 

The  true  shape  of  the  beam  envelope  and  of  the  ion  streamlines  can  be 
revealed  by  inspecting  a set  of  radial  density  profiles  scanned  at  various  axial 
positions.  Such  radial  scans  have  been  obtained  using  cylindrical  probes  as  shown 
on  Fig.  III-4,  and  also  using  small  collimated  Faraday  cups  as  reproduced  from 
experimental  traces  in  Fig.  III-5.  Only  marginal  misalignments  of  the  source  with 
respect  to  the  tunnel  centre-line  are  observed,  as  indicated  in  Fig.  III-4,  and 
are  typically  less  than  1.5°  • The  slight  distortions  which  could  have  been  intro- 
duced on  the  radial  profiles  by  keeping  the  probes  parallel  to  the  centre-line, 
rather  than  following  an  arc  to  preserve  the  true  distance  from  the  source,  are 
not  believed  to  be  significant  because  of  the  small  corrections  involved.  The  small 
envelope  irregularities  are  more  easily  observed  with  the  Faraday  cups  as  shown  in 
Fig.  III-5  and  are  more  pronounced  at  the  extreme  edge  of  the  beam.*  One  can  use 
these  profiles  to  obtain  an  idea  of  the  departure  of  the  ion  streamlines  from  the 
straight  lines  assumed  by  the  conical  model.  This  is  accomplished  by  taking  the 
width  of  the  profile  at  one  e-fold  from  its  maximum  (~37 %)  since  the  beam  structure 
is  found  to  have  the  familiar  bell- shape  feature  of  a Gaussian  model.  'The  widths 
are  obtained  directly  from  the  normalized  ion  density  profiles  which  are  measured 
experimentally,  as  depicted  in  Fig.  III-6.  The  various  radial  positions  are 
normalized  with  the  corresponding  ideal  source-flow  values  of  the  beam  radius  rB 
at  an  axial  position  z,  that  is 

r = r + ztan7  (ill.  21) 

Jd  O 

The  half-width-at-one-e-fold  (HWOE)  measurements  can  be  plotted  against  the  non- 
dimensional  position  z/rQ  for  a number  of  beam  conditions  as  shown  in  Fig.  III-7. 

The  results  clearly  indicate  the  gradual  bending  of  the  ion  streamlines  up  to 
axial  distances  of  the  order  of  30  for  the  large  beam  diameter,  and  of  roughly 
z/rD  ~ 60  in  the  case  of  the  small  beam  diameter . The  expansions  become  basically 
straight  or  conical  beyond  these  points  and  follow  closely  the  measured  spreading 
angles  7 in  each  case.  The  fact  that  the  various  beam  conditions  do  not  lead  to 
a unique  result, as  would  be  expected,  required  some  explanation.  This  fact  was 
already  observed  in  the  normalized  profiles  presented  in  Fig.  II-6  where  the  two 
solid  curves  represent  an  attempt  to  average  the  three  z/r0  cases  for  both  the 
small  and  the  large  beam  diameter  configurations.  The  answer  to  this  observation 
rests  in  the  experimental  traces  themselves.  The  small  beam  exit  diameter  will 
generally  lead  to  relatively  flatter  profiles  because  the  beam  feeds  from  a smaller, 

*These  anomalies  could  be  responsible  for  a somewhat  irregular  probe  behaviour 
reported  in  Section  17.4.2. 
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more  uniform  portion  of  the  discharge.  The  plasma  density  is  also  considerably 
smaller  and  it  is  therefore  more  difficult  to  extract  the  signal  from  the  traces, 
as  indicated  by  the  larger  spread  of  the  solid  symbols  in  Fig.  III-6,  especially 
for  the  larger  axial  distances  where  the  signal  level  is  lowest.  Furthermore, 
a small  background  contribution  of  slow  ions  has  not  been  subtracted  from  these 
traces;  it  constitutes  a relatively  larger  contribution  to  the  lower  signal. 
Finally,  a further  contribution  to  the  deviation  can  be  due  to  the  fact  that  the 
density  profiles  are  not  perfectly  similar  in  shape.  Here  we  have  treated  them 
as  ideal  Gaussians,  according  to  the  standard  representation 


n(r)  = n(o)  exp 


(III. 22) 


A final  check  of  the  beam  density  model  is  provided  by  a series  of  axial 
scans  obtained  with  a group  of  cylindrical  probes  perpendicular  to  the  flow  at 
roughly  the  same  beam  energy.  The  density  was  measured  using  the  perpendicular 
probe  relation  obtained  in  Chapter  I 


J 


en.  • d £ • 7/  1 + <X>  ' 
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(III. 23) 


where  $p  is  the  nondimensional  probe  bias  and  after  the  contribution  of  the  slow- 
ions  is  subtracted  from  the  measured  current  according  to  a method  outlined  in 
Section  IV .1.1.  The  results  are  presented  in  Fig.  III-8  where  ni  is  plotted 
against  z/rQ  and  is  normalized  with  respect  to  the  value  at  z/rQ  = 20  and  4o  for 
the  large  and  the  small  beam  diameter  cases,  respectively.  The  conical  model 
lines  which  satisfy  the  ideal  source-flow  relation 


ni 


(111.24) 


are  drawn  to  provide  a direct  comparison.  The  trend  indicated  corroborates  the 
previous  experimental  observations,  namely  that  the  ion  streamlines  bend  slightly 
outward.  This  suggests  the  existence  of  an  axial  beam  density  fall-off  which 
should  be  slightly  faster  than  (z/rG)"2  at  the  larger  distances,  as  indicated  in 
Fig.  III-8 . 

In  summary  then,  we  have  observed  that  the  ideal  source-flow  model  is 
somewhat  distorted  due  to  the  slight  curvature  of  the  ion  streamlines  within  the 
beam  envelope.  This  is  also  reflected  by  an  ion  density  which  decreases  slightly 
faster  than  (z/rQ)  2 at  large  axial  distances.  The  immediate  consequences  for 
the  geometrical  model  used  to  describe  the  simulation  of  the  transverse  teuperature 
are  difficult  to  delineate  quantitatively.  However,  since  Si  is  intimately 
related  to  the  assumed  model  for  the  beam  formation  and  ®eo  depends  on  7, 

u>  z/ro  an<i  9 according  to  Eq.  III. 13  and  Eq.  III-16,  we  can  expect  that  the  speed 
ratio  simulated  will  be_ somewhat  perturbed  by  any  significant  modification  of  the 
actual  values  of  7 and  u.  For  instance,  the  value  of  the  divergence  angle  7 is 
decreased  at  small  z/rQ  because  of  the  existence  of  a beam  waist.  Therefore  the 
limiting  transverse  velocity  defined  by  III. 13  could  be  decreased  at  some  z/rOJ 
thus  the  synthetic  speed  ratio  would  consequently  increase.  Finally,  at  large 
z/r0,  larger  transverse  velocities  are  expected  because  the  ion  streamlines  are 
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curved  outward.  This  effect  would  lead  to  speed  ratios  slightly  overestimated 
by  the  ideal  conical  expansion,  as  discussed  in  Chapter  IV,  Section  IV. 3.1. 


Ill .2  Pertinent  Plasma  Features 
III. 2.1  Background  Ions 

The  high-energy  ions  constitute  a fraction  g of  the  total  gas  flux  from 
the  accelerator.  The  remaining  portion  comes  from  the  non-ionized  neutral  parent- 
gas  such  that 


g*T—-r+-r  (III.25) 

tot  + a . 

where  the  neutral  flux  is  roughly  given  by  the  one- dimensional  Maxwellian  gas 
relation 


and 


n „ anode 

raa 


(III  .2 6) 


(III. 27) 


A background  of  slow  ions  is  produced  by  the  charge- transfer  mechanism  schematically 
represented  by  the  relation 


A+  + A A+  + A (111.28) 


This  mechanism  is  responsible,  in  the  low-density  plasma  facilities  of  large 
dimensions,  for  the  creation  of  an  annoying  fraction  of  slow  randomized  ions  whereby 
an  energetic  ion  loses  its  charge  to  a slow  neutral  atom,  indicated  by  the  symbol 
~,  following  a rate  process  governed  by  a cross-section  QcT  such  that 


n. 

,n7  = 1 ‘ ex^-  ^cT  na  ~ ^cT  na  2 (XII.29) 

The  axial  distance  travelled  z can  be  significant  for  extended  test-sections.  Thus, 
many  authors  (Refs.  19,24,30)  have  reported  relative  values  for  slow-ion  densities 

n. 

0 < — < 20$ 
n. 


for  facilities  using  an  ion  accelerator  similar  to  the  one  used  here.  Using  a 
method  obtained  from  the  work  of  Clayden  and  Hurdle  (Ref.  29)  we  have  monitored 


r-”  " — - - 

the  slow-ion  contribution  within  the  beam  using  small  flat  cups  (Section  II. 4). 
These  measurements  have  shown  that  the  maximum  relative  slow-ion  density  within 
the  beam  favours  the  higher  figure  reported  above,  when  the  small  cups  were 
i looking  back  on  axis  within  the  beam,  whereas  we  get 


n. 

1 


with  the  cups  looking  across  the  beam.  The  measured  values  also  indicate  that 
the  relative  background  density  scales  roughly  as 


and 


which  agrees  with  the  fiist  order  estimate  obtained  from  Eqs.  III. 23,  25,  27  and 
the  conical  density  relation  given  in  111.24,  i.e. 


/v 


(III. 30) 


since  QcT  is  not  varying  much  (Ref.  72). 

In  order  to  evaluate  the  importance  of  the  slow-iors  which  introduce  an 
offset  current  independent  of  probe  orientation  on  our  probe  traces  (Fig.  IV-l) 
we  have  obtained  a background  pressure  dependence  by  introducing  a controlled  back- 
fill of  neutral  argon  in  the  test  chamber.  This  was  judged  preferable  rather  than 
feeding  the  gas  through  the  ion  source,  as  it  should  be,  so  that  the  operating 
conditions  of  the  accelerator  would  not  be  altered.  Typical  results  are  presented 
in  Fig.  III-9  where  the  background  flow  ion  current  is  displayed  in  arbitrary  units. 
If  we  assume  quite  reasonably  that  only  a marginal  partial  pressure  of  foreign  gases 
is  present  in  the  test  section  during  normal  operation,  we  can  see  that,  according 
to  the  ideal  gas  law  and  to  Eq.  III. 29  we  should  have 


n. 

-oCn 
n.  a 
i 


of  P 


where  P is  the  background  pressure.  This  is_verified  in  Fig.  III-9  except  for 
values  of  test  section  pressures  near  8 x 10  5 torr  beyond  which  back- ingest ion 
into  the  source  produces  severe  beam  instability.  By  directly  extrapolating  this 
linear  region  we  estimate  that  the  residual  slow-ion  contribution  would  become 
negligible  at  about  8 x 10  6 torr  which  is  somewhat  lower  than  our  minimum  operating 
point . 
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The  background  ion  current  is  superposed  on  any  measurement  obtained  by 
the  cylindrical  probes.  This  very  fact  has  led  us  to  use  another  single  method  to 
estimate  the  operating  pressure  at  which  the  stray  contribution  of  randomized  slow 
ions  would  be  negligible.  This  is  done  by  looking  at  the  current  collected  by  the 
perpendicular  probes  (0  = 7r/2) . The  absolute  current  to  the  probe  in  the  ideal 
situation,  that  is  without  any  slow  ions,  is  given  by  the  following  relation 
obtained  from  1.7 


J - , <t> 

80  \ 2 5 p 
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(1  + 4>  )2 
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where 
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= en.  u 
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Therefore,  by  plotting  the  ratio  of  the  measured  value  and  the  one  given  above 
versus  background  pressure,  one  obtains  an  estimate  of  the  spurious  contribution^ 
A group  of  5 experimental  cases  are  presented  in  Fig.  III-10  for  2 probes,  at 
about  the  same  beam  energy  value  and  for  three  axial  positions.  The  pressure 
exclusion  was  achieved  as  before  using  a backfill  of  neutral  argon.  A single 
extrapolation  yields  about  1 x 10_5  torr  for  the  pressure  for  negligible  slow  ion 
contributions  which  is  quite  close  to  the  pressure  estimated  with  III. 30.  By 
inspection  we  can  then  verify  that  the  relative  importance  of  the  stray  contribu- 
tion does  not  exceed  20$  at  the  typical  operating  conditions. 


III.2.2  Plasma  Potential 

Clayden  (Ref.  29)  has  obtained  the  dynamic  equivalent  to  Chen's  'static' 
relation  for  the  estimate  of  the  floating  potential  in  a plasma  (Ref.  73) • These 
relations  yield  basically  identical  results  for  a dynamic  plasma  with  Te  ~ l/3  eV 
and  with  a directed  ion  energy  of  500  eV,  that  is 


( Chen- static) 
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(Clayden) 


where 
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(III.32) 


This  has  led  us  to  consider  the  value  of  Vf  measured  with  respect  to  reference 
ground  as  an  approximated  value  of  the  local  plasma  potential,  since  Vf.  is 
typically  close  to  50  volts.  Therefore,  with  an  accuracy  of  about  2$,  one  can 
obtain  the  nondimen si onal  probe  bias 

0.  - 0„ 


0 ■ 

P \ Mu2 


31 


by  using  the  following  approximations 


and 


\ M^2  = eV.  = e(VB  - e(V+  + VD  - Vf) 


We  measure  Vf  by  taking  the  reading  of  the  bias  at  which  the  probe  collects  no 
current  when  placed  within  the  plasma.  The  characteristic  values  of  the  floating 
potential  are  shown  in  Fig.  III-ll  as  a function  of  beam  energy,  for  number  of 
probes.  The  relatively  large  values  of  floating  potential  (with  respect  to 
reference  ground)  which  are  achieved  are  typical  for  the  operation  of  ion  acceler- 
ators of  this  type  when  the  neutralizer  wire  is  withdrawn  from  the  beam  (Ref.  65). 
In  fact,  in  these  experiments  the  neutralizer  was  placed  outside  the  nominal  beam 
diameter.  In  the  low- voltage  configuration,  indicated  by  solid  circles  in  the 
figure,  the  beam  fans  out  prohibitively  at  the  lower  beam- forming  voltages.  This 
has  basically  the  same  effect  as  that  of  re-inserting  the  neutralizer  within  the 
beam  thereby  decreasing  the  coupling  voltage  between  it  (or  reference  ground) 
and  the  plasma. 


\ 
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III. 2. 3 Electron  Temperature 

Although  the  absolute  value  of  the  electron  temperature  is  not  used 
specifically  in  the  analysis  of  the  probe  response  which  is  pertinent  to  this 
report,  an  order  of  magnitude  estimate  of  that  value  was  judged  essential  to 
have  a 'complete'  picture  of  the  plasma.  Three  different  probe  geometries  were 
used.  Basically  all  three  methods  were  applied  while  assuming  static  conditions. 
The  standard  procedures  obtained  from  the  literature  were  followed  for  the  single- 
probe, the  loop-probe  and  for  the  double-probe  - namely,  using  the  slope  of  a 
semi-log  plot  of  current  versus  voltage  (Refs.  5,36,73)-  An  approximated  tangent 
method  described  in  Appendix  E was  also  used.  Although  no  attempt  is  made  in 
order  to  favour  one  method  over  another,  one  can  say  that  the  loop  has  the 
distinctive  advantage  of  being  the  only  self-cleaning  (ohmic)  candidate  of 
interest.  A table  shown  in  Fig.  Ill- 12  summarizes  the  various  results  which  are 
presented  in  Fig.  III-13-  Only  a marginal  dependence  on  axial  position  is 
detected  by  all  probes.  The  results  are  basically  independent  of  probe  position 
(perpendicular  or  parallel) . The  loop-probe  results  are  slightly  higher  than  the 
neutralizer  wire  temperature 

Tg  (loop)  ~ 0.35  eV 
Tw  ~ 0.27  eV 

The  loop  geometry  is  judged  more  compatible  with  the  'static'  plasma  approach  as 
the  cylindrical  single-probe  features  are  more  liable  to  be  distorted  by  dynamic 
effects,  such  as  the  end  effect  outlined  in  this  report  (see  Appendix  J) . The 


J 
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results  obtained  with  the  single-probe  axe  systematically  about  100$  higher  than 
Te  (loop).  The  approximation  scheme  mentioned  in  Fig.  III-13  refers  to  a descrip- 
tion reported  in  Appendix  K.  At  any  rate,  a value  of  about  l/3  eV  is  taken  as 
representative  of  the  electron  temperature  for  the  standard  beam  conditions  reported 
here . 


EXPERIMENTAL  RESULTS  AND  DISCUSSIONS 


IV. 1 Experimental  Procedures 


IV. 1.1 


jerimental  Probe  Response 


After  regular  and  systematic  probe  cleaning  (see  Ch.  II. 3 , Appendix  C) 
the  cylindrical  probes  were  positioned  by  adjusting  the  X and  Y coordinates  until 
a peak  reading  of  the  ion  current  was  achieved  with  the  probe  parallel  to  the  flow 
geometric  centre-line.  The  sensitive  pitch  and  yaw  adjustment  actuators  were  then 
used  to  further  optimize  the  peak.  This  procedure  effectively  yielded  the  local 
direction  of  the  ion  velocity  because  of  the  probe  attitude  sensitivity,  and  at 
this  point  the  recording  of  the  various  flow  and  probe  parameters  was  initiated. 


The  salient  features  of  the  ion  current  versus  pitch-angle  curve  are 
exemplified  in  Fig.  IV-1.  The  solid  line  represents  the  actual  finite  probe  current 
response  for  values  of  0 from  0°  to  90°.  As  the  theory  suggests,  the  ion  current 
collection  is  expected  to  depart  significantly  from  the  'infinite'  probe  behaviour 
for  small  angles  only.  The  dotted  line  thus  represents  the  current  extrapolated 
at  9 = 0°  which  this  length  of  probe  would  have  collected,  under  the  same  conditions, 
if  it  would  have  been  part  of  an  infinitely  long  cylinder.  This  extrapolation  is 
achieved  by  fitting  a flexible  curve  tracer  or  spline  to  the  curve  for  pitch  angles 
between  about  ± 45°  about  the  peak  and  by  smoothly  joining  the  two  sides  of  the 
narrow  peak.  The  values  of  the  'infinite'  probe  current  at  0 = 0,  .1^(0)  estimated 
this  way  were  found  to  be  very  close  to  the  calculated  values  from  Eq.  1.4  where 
the  current  at  0 = 90°,  J (tt/2)  is  basically  the  'infinite'  probe  result.  This 
simple  method  was  therefore  used  with  confidence  on  all  the  experimental  traces. 

Since  the  current  peak  is  determined  by  the  fast  beam-ions  only,  the  true  base- 
line for  the  given  experimental  conditions  shown  dotted  here  is  obtained  when  the 
slow-ion  background  contribution  J is  subtracted  uniformly  for  all  angles,  as 
discussed  later.  The  full-width-at-half-maximum  (FWHM)  0^  is  obtained  by  measuring 
the  angular  width  of  the  peak  at  half  its  height  measured2f rom  the  true  base-line. 
Another  half-width  01,  defined  this  time  by  the  full  width  at  the  half  distance  between 
the  peak  maximum  J(0^  and  the  'infinite'  probe  extrapolation  (0),  was  used  to  com- 
pare the  results  with  the  analysis  of  Hester  and  Sonin  (Ref.  18) . Typical  traces 
for  four  values  of  probe  voltage  are  shown  on  Fig.  1-3  together  with  an  angle  template 
which  was  drawn  as  an  overlay  for  convenience. 


The  basic  flow  conditions  together  with  the  probe  response  features 
described  above  were  measured  systematically  on  the  experimental  traces  and  re- 
corded on  perforated  cards  for  data  reduction  using  a computer  code. 

As  was  explained  earlier  in  Chapter  III. 2. 2,  the  local  plasma  potential 
can  be  estimated  by  using  the  floating  (zero-current)  potential  Vf  as  shown  on 
Fig.  IV- 2 to  yield  the  nondimensional  probe  bias  with  respect  to  system  ground 

e(V  - V ) V - V„ 

<D  = . P space  „ jo f 

P i " VB  - Vf 


, ....  . . 


where  V is  the  effective  accelerator  voltage 

D 


V = V,  + VT 


where  Vp  is  the  discharge  potential  and  V+  and  V_  the  voltage  on  the  extracting 
and  the  accelerating  grid,  respectively,  A slight  difference  of  about  4 volts 
is  typically  found  between  Vf  for  the  perpendicular  and  the  parallel  probe 
orientations,  and  an  average  value  is  satisfactory  in  view  of  the  relatively 
large  probe  biases. 


The  background  (slow-ion)  ion  current  J can  be  estimated  and  subtracted 
from  the  trace  as  follows.  We  make  the  reasonable  assumption  that  the  slow-ion 
spurious  current  is  collected  by  the  probe  without  any  attitude  sensitivity,  as 
if  it  were  in  a static  plasma.  Thus,  using 


we  get 


jjf)  _ 

Jjoj 


i + $ 


= K (no  slow-ions) 
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1 ocA  2 ^ 


extrapolated 

(J  (£)  + J)  " (J  (0)  + J)  ' K i 
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where  J (0)  is  the  'infinite-probe'  current  at  0 = 0 obtained  by  extrapolation. 
Each  experimental  probe  record  is  treated  in  this  systematic  way  to  acquire  the 
peak  features  and  the  flow- field  characteristics.  To  minimize  the  effect  of 
varying  parameters  not  directly  pertinent  to  the  experiment,  a set  of  standard 
source  conditions  were  rigorously  followed  through  the  extent  of  this  study. 


IV. 1.2  Range  of  Parameters 


As  the  theory  has  suggested,  the  characteristics  of  the  end  effect 
depend  on  a significant  number  of  parameters  which  have  to  be  accounted  for  to 
be  able  to  measure  the  ion  speed  ratio  or  the  ion  temperature  using  this  effect. 
The  most  pertinent  parameters  of  interest  are  those  of  the  plasma  stream  itself, 
namely  the  ion  density,  the  ion  bulk  velocity,  the  ion  speed  ratio  and  the 
electron  temperature  and  those  related  to  the  presence  of  the  probe,  that  is, 
the  probe  bias,  the  probe  aspect  ratio  ($/ d)  and  the  nondimens ional  probe  length 
Tg  which  is  the  major  probe  parameter  when  characterizing  the  finite-length 
effect . 


Following  our  discussions  in  Section  III  we  can  attempt  to  tailor  the 
ion  speed  ratio  by  using  a simple  model  developed  for  a conical  beam  geometry. 
The  concept  of  synthetic  temperature  introduced  in  Section  III  suggests  that  the 
ion  speed  ratio  would  scale  directly  with  the  nondimensionalized  axial  distance 
z/rQ  and  would  depend  only  on  z/rD.  Therefore  both  the  axial  distance  z and 
the  exit  beam  radius  rQ  were  varied  to  cover  a range  of  speed  ratios  from  about 


10  to  100.  Ion  beam  directed  energies  between  60  eV  and  1200  eV  were  considered. 
Since  the  end  effect  is  not  expected  to  have  any  first  order  sensitivity  to  the 
electron  temperature,  no  attempt  was  made  to  vary  this  parameter  within  the 
plasma  flow. 


The  probe  parameters  discussed  in  our  theoretical  considerations 
(Chapter  I)  are  basically  dictated  by  the  range  of  validity  of  the  analysis  for 
the  end  effect.  In  order  to  cover  the  widest  range  possible  for  the  probe  bias, 
was  varied  from  about  0.03  to  1.  For  the  nominal  operating  regime  of  the 
plasma  source  the  conical  beam  ion  density  varies  roughly  as  (z;/r0)~2.  Probe 
diameters  of  0.01  to  0.05  cm  were  used  to  test  the  sensitivity  of  the  effect 
to  the  probe  radius  to  Debye  length  ratio,  6,  within  a range  of  about  0.01  to 
0.5.  Finally,  since  the  finite  length  effect  is  most  noticeable  for  values  of 

Ti’ 


smaller  than( about  1.5  (Ch.  I),  this  requirement  dictated  the  use  of  probe  £/d 
ratios  from  about  20-500,  such  that  we  eventually  covered  To  values  from  0.01 
to  1.5.  Figure  II- 11(b)  gives  all  the  experimental  runs  according  to  the 
various  ranges  outlined  above. 


IV .2  Verification  of  the  Finite-Probe  Theory 

IV. 2.1  .Behaviour  of  the  Main  Peak  Features 

A first  basic  inspection  of  Sanmartin’s  theory  can  be  made  by  the 
determination  of  the  dependence  of  the  probe  current  response  on  probe  bias. 

The  two  main  features  of  interest  concerning  the  peak  are  undoubtedly  its  height 
and  its  width  since  they  both  convey  valuable  ion  temperature  information.  The 
analysis  indicates  that  the  peak  height  should  grow,  to  first  order,  as 

j(e  =0)  ar$p  (rv.i) 

and  that  the  width  at  half-maximum  should  basically  be  independent  of  <£p. 

A typical  example  of  this  dependence  is  shown  on  Fig.  IV. 3 where  a 
series  of  traces  are  superposed  to  show  the  peak  growth  with  varying  probe  bias. 
The  peak  is  vanishingly  small  for  values  of  probe  voltage  near  floating  potential. 
A slight  asymmetry  which  can  be  seen  on  the  extrapolation  at  0 = 0 is  ascribed 
to  not  quite  ideal  alignment  of  the  probe  on  the  centre-line.  Once  the  slow-ion 
contribution  is  properly  subtracted,  the  linear  growth  of  the  peak  value  with 
absolute  probe  bias  can  be  verified,  as  shown  on  Fig.  IV-4  for  a number  of  probes. 
One  notices,  in  particular,  the  extent  of  probe  bias  covered  which  indicates  that 
the  linear  relationship  is  preserved  for  very  large  voltages.  The  peak  half- 
width was  found  to  be  independent  of  the  bias  in  accordance  with  the  theory,  for 
nondimensional  probe  bias  values  as  high  as  l/2  as  shown  on  Fig.  IV-5  for  6 
different  experiments.  In  addition,  the  product 

7-  (9  = 0)  • /Qi  of  (IV.2) 
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is  found  to  scale  according  to  the  theory,  confirming  the  proportionality  of 
J (0  =0)  with 

00'  ' p 

IV .2 .2  Comparison  With  Previous  Experiments 

The  dependence  of  the  current  to  the  finite  parallel  cylindrical 
probe  on  the  nondimensional  probe  length  rg,  as  first  proposed  by  Hester  and 
Sonin  (Ref.  18) , was  corroborated  by  these  authors  themselves  with  a series  of 
experiments  conducted  with  three  probes  in  an  argon  beam.  The  ion  accelerator 
they  used  was  an  electron-bombardment  ion  thruster  very  similar  to  the  one 
used  here.  Hester  and  Sonin  report  the  following  experimental  conditions: 

d = 0.013  cm,  ^ = 75,135,225 

^p  < "15 

# 

5 < Ma  < 40 

T. 

4000°K  < T < 60,000°K  and  ^ « 1 
e 5 T 

e 

0.2  < T,  < 4 

As  a matter  of  comparison  these  results  and  the  analytical  curves 
pertinent  to  the  range  of  parameters  studied  by  Hester  and  Sonin  are  reproduced 
from  Sanmartin  (Ref.  22)  on  Fig.  IV-6,  as  a function  of  probe  length  and  the 
ratio  of  probe  radius  to  Debye  length.  Two  curves  obtained  by  Sanmartin  for 
small  tenqperature  ratios  Tj_/Te  are  shown  along  with  the  curves  of  Hester  and 
Sonin  which  are  for  zero  values  of  this  ratio.  The  results  agree  relatively 
well  with  the  theoretical  predictions. 

In  order  to  produce  some  comparison  with  the  Hester  and  Sonin  numerical 
solutions  (Ref.  18)  we  have  shown  on  Fig.  IV-7  two  groups  of  experimental  results 
obtained  in  the  current  investigation  for  a range  of  temperature  ratio  signif- 
icantly larger  than  the  ones  reported  previously, 


0.5  < P < 3*5  (this  report) 

B < 2 x 10  2 (Hester-Sonin) 


together  with  3 theoretical  curves  obtained  from  Hester  and  Sonin  (Ref.  18) . The 
experiments  were  sorted  into  two  groups,  which  are  distinguished  by  different 
ranges  for  the  p ratios.  Since  the  experimental  results  within  each  group  were 
obtained  for  a relatively  wide  range  of  experimental  conditions,  they  are  not 
expected  to  belong  to  a single  P-e  curve,  however  we  wish  to  indicate  that  the 
Hester-Sonin  predictions  are  relatively  well  compatible  with  the  experimental 
results  obtained  for  the  high  p values  reported  here. 
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Both  the  numerical  computations  provided  by  Hester  and  Sonin  and 
Sanmartin's  more  elaborate  analytical  description  demonstrate  that  the  end 
effect  would  not  display  any  first  order  sensitivity  to  electron  temperature. 

A graphical  representation  such  as  the  one  on  Figs.  IV- 6 and  IV-7  however 
may  seem  a bit  confusing  to  the  reader.  The  experimental  curves  shown  dotted 
on  Fig.  IV-6  for  the  two  values  of  temperature  ratio  (£)  of  10  2 and  10  3 appear 
to  demonstrate  a sensitivity  to  electron  temperature  which  could  be  significant 
for  the'  smaller  values  of  e,  that  is,  for  the  range  where  the  sheath  radius  or 
the  Debye  length  are  much  larger  than  the  probe  radius.  In  fact,  we  know  that 
Hester  and  Sonin  have  varied  the  absolute  value  of  Te  from  about  0.3  to  5 eV. 

The  apparent  effects  of  P and  e are  deceptive  since  a change  in  Te  varies  the 
values  of  P,  e and  ip^  simultaneously  in  such  a way  as  to  keep  j/ JodtCo)  nearly 
constant.  Therefore  we  feel  that  the  graphical  representation  adopted  is  not 
ideal.  The  contention  that  the  ion  current  peak:  structure  is  not  influenced 
significantly  by  the  electron  temperature  is  corroborated  more  clearly  by  the 
more  appropriate  choice  of  parameters  we  have  introduced  in  Chapter  I.  The 
dependence  on  Te  cancels  out  except  by  the  way  of  6 which  only  involves  a weak 
logarithmic  dependence  (Appendix  K) . One  can,  in  fact,  account  for  the  correct 
description  of  the  end  effect  by  using  a set  of  nondimensionalizing  parameters 
which  do  not  depend  on  T_,  as  we  will  see  in  the  next  sections  (Figs.  IV-10, 

IV- 11) . 


TV. 2. 3 Self-Consistent  Measurement  of  the  Speed  Ratio 

The  method  we  have  derived  from  Sanmartin's  analysis  to  determine  the 
ion  speed  ratio  using  the  end  effect  stems  from  the  speed  sensitivity  of  both 
the  peak  height  and  the  peak  width.  This  means  that  a measurement  of  the  current 
peak  height  at  0 = 0 and  of  the  angular  width  at  half  the  peak  height  can 
reinforce  one  another  in  determining  the  ion  speed  ratio.  This  is  accomplished 
by  using  both  measurements  obtained  from  a single  trace  to  determine  the  value 
of  the  parameter  5,  which  is  the  key  to  obtain  S^  as  was  indicated  in  Section 
1.1.3-  The  interdependence  of  peak  height  and  width  with  the  various  parameters 
and  a schematic  representation  of  the  interpretation  procedure  is  given  in  Fig. 
IV-8.  A value  of  5 is  determined  by  iteration  locating  the  intersection  of  two 
functions,  as  demonstrated  on  Fig.  IV-9  for  two  experimental  cases. 


The  procedure  starts  with  the  measured  quantities  (which  dete:  line  the 
value  of  the  peak- shape  function  X obtained  by  two  explicit  relations  involving 
the  unknown  value  of  5: 


X = X(J (0),  5,  4>p,R  , x g)  = X jl  (5) 

X = X(0i,  6,  0 R , t.)  = X 2 (6) 

2 P * 


( iv. 3) 


Therefore,  there  must  exist  a 8 hereafter  called  the  self-consistent  5 such  that 


Xl^selP  Xs(Jself) 
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Upon  finding  the  solution,  both  6 and  X are  known  and ‘the  value  of  Si  is  obtained 
from  the  analytical  expressions  (Ch.  I,  Eqs.  1.15  to  I.l8).  Since  X(0  = 0) 
depends  on  the  parameter  t) , this  latter  parameter  can  in  turn  be  determined  from 
Eq.  1.19  which  has  the  general  form 


f 

fc 


= ^6self’  V R’  Ti*  Si} 

and  has  to  be  solved  numerically.  A simple  iterative  computer  program  was  used 
to  determine  Si  using  the  set  of  data  points  obtained  from  all  the  recorded  traces 
using  this  scheme. 


.IV. 2. 4 Sensitivity  of  the  Peak  Features  to  Speed  Ratio 

In  order  to  demonstrate  the  ability  to  describe  the  end  effect  completely' 
using  the  proposed  set  of  parameters  which  are  readily  available  from  measurement, 


V R>  Tr  V 5 

we  have  presented  on  Figs.  IV-10  and  IV- 11  the  solution  for  the  peak  height  and 
the  peak,  width  in  terms  of  these  parameters.  As  a global  indication  of  experimental 
behaviour,  data  points  for  a large  number  of  experiments  are  shown  sorted  into  four 
speed  ratio  ranges.  These  experiments  all  correspond  to  about  the  same  value  of 
the  parameter  e*.  This  is  not  a parameter  containing  extra  information  as  we 
have  seen,  since  it  represents  the  ratio  of  the  two  nond.imensional  probe  lengths 
R and  Tg 


e*  s 


1L. 

n/2  R 


This  ratio  in  fact  conveys  the  information  about  the  relative  size  of  the  probe 
radius  in  terms  of  the  length- scale  defined  in  a fashion  akin  to  the  electron 
Debye  length  but  containing  the  directed  ion  energy  instead  of  the  electron 
thermal  energy,  i.e. 


e* 


r 

= 


and  oQ  = 


The  theoretical  curves  in  Figs.  IV-10  and  IV-11  are  drawn  for  convenience  with  the 
following  set  of  values 


e*  = 0.001 


0 =0.10 

P 

6 = 0.25 
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and  the  nondimensional  peak  current  results  are  scaled  according  to  the  relation- 
ship we  have  verified  previously 


(0  = 0)  of  6 

J ' P 

00  * 

and  plotted  with  the  theory  whereas  the  peak  width  results  also  obtained  for  the 
range  of  probe  bias 


0.03  < <t>  <0.5 
P 

were  plotted  directly  since  they  we  re  found  to  be  independent  of  both  6 and  $ . 
The  four  groups  of  data  points  include  measured  speed  ratios  near  five  nominal! 
values,  namely 


S±  = 20,  30,  4o,  60,  80 

obtained  for  various  probe  lengths  and  diameters.  The  agreement  is  excellent  for 
0.03  < Tg  < 1.  The  experimental  range  of  Tg  was  selected  for  values  smaller  than 
about  1.5  for  reasons  evident  in  Figs.  IV-10  and  TV-11.  In  fact,  the  peak  width 
levels  off  for  values  of  smaller  than  1,  in  general,  whereas  the  peak  height 
reaches  a maximum  at  about  1 and  decreases  slowly  thereafter.  This  region  of 
interest  is  fortunately  well  suited  for  the  range  of  operation  of  the  ion  source 
we  used. 


Finally,  in  order  to  conclude  this  verification  of  Sanmartan' s theory 
concerning  the  dependence  of  the  peak  features  on  and  we  can  compare  the 
peak  half-width  measured  from  the  traces  with  the  half- width  predicted  by  Bettinger 
and  Chen's  theory  (Ref.  17).  Using  the  measurements  of  Te  (Fig.  III-13)  and  the 
estimates  of  the  ion  temperature  we  are  simulating, using  the  synthetic  temperature 
model  developed  in  Chapter  III,  che  temperature  ratio  T^/Tg  pertinent  to  this 
investigation,  is  found  to  be  of  order  of  unity.  The  value  of  01  obtained  from 
the  end  effect  theory  of  Bettinger  and  Chen  can  therefore  be  evaluated: 

T. 

At  p = — 
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(IV. 5) 


An  estimate  of  the  trend  can  be  obtained  from  Fig.  IV-12  on  which  we  have  plotted 
the  ratio  of  the  experimentally  measured  half- width  and  the  half- width  from  Eq. 
IV.  5: 

/v 

0i_  (measured)/©^ 

2 2BC 


39 


We  have  added  the  dotted  lines  for  three  prohe  i/d  (r)  values  as  a matter  of 
comparison  using  the  analytical  expression  of  91.  obtained  from  Sanmartin  with 
<t>p  =0.1  and  5 = 0.25.  The  overall  behaviour  indicates  that  in  this  range  of 
values  for  t£  the  peak  half- widths  are  not  correctly  predicted  by  Bettinger  and 
Chen’s  analysis.  At  large  rg  the  ratio  levels  off  to  values  significantly  larger 
than  unity,  in  contrast  with  the  results  obtained  for  ion  to  electron  temperature 
ratios  near  zero,  reported  previously  (Ref.  l8) . This  result  for  T^/Tg  values 
near  unity  is  consistent  with  the  theoretical  values  predicted  by  Sanmartin,  as 
was  the  one  reported  by  Hester  and  Sonin,  and  shows  that  for  an  ion  speed  ratio 
(Si)  of  10  with  a typical  e*  ~ 0.005,  then  for  large  Tg  and  T^/Tg  ~ 1, 
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IV. 3 Discussion  of  Experimental  Results 

IV.3.1  Verification  of  the  Geometrical  Speed  Ratio  Model 

In  order  to  assess  the  merit  of  the  method  developed  to  measure  the 
ion  speed  ratio  using  the  end  effect  of  cylindrical  probes,  we  have  to  confront 
the  results  with  the  model  predictions.  A group  of  experimental  results  repre- 
sentative of  the  range  of  parameters  studied  are  presented  on  Fig.  IV-13  to 
exemplify  the  comparison  between  the  self-consistent  values  of  with  the 
model-line  which  is  drawn  according  to  the  simple  relation 


o =3  tt 

si  ~B~ 

geo 


The  trend  seems  to  favour  a power  law  which  is  slower  than  (z/rD)+  . The  overall 
behaviour  therefore  indicates  that  the  measured  speed  ratios  are  higher  than  the 
model  suggests,  whereas  at  the  high  z/rQ  the  values  of  Sj_  measured  fall  short  of 
the  model.  The  results  obtained  with  the  very  short  probe  (open  diamonds)  suggest 
that  the  probe  behaviour  for  small  aspect  ratios  differs  appreciably  from  the 
theoretical  predictions.  Specific  conclusions  in  this  regard  are  difficult  because 
the  sensitivity  of  such  short  probes  to  speed  ratios  is  predicted  to  be  low  for 
values  of  above  about  30  as  shown  in  Figs.  1-7  and  1-8. 


This  general  departure  from  the  geometrical  model  can  be  explained  in 
terms  of  real  flow-field  effects  rather  than  probe  effects.  The  arguments  which 
justify  the  synthetic  ion  temperature  simulation  are  based  on  the  simplified 
geometry  of  an  ideal  conical  beam.  In  fact,  both  the  radial  scans  and  the  photo- 
graphic records  have  led  to  observations  indicating  that  the  beam  was  somewhat 
distorted  from  the  conical  model  (see  Figs.  11-10  and  III-10)  as  is  fairly  typical 
of  ion  streams  emitted  from  small- size  ion  thrusters  (see  Refs.  27-30  on  the 
subject).  A complex  combination  of  plasma-electrostatic  field  couplings  in  the 
vicinity  of  the  beam  shaping  grids  is  responsible  for  the  definition  and  the 
stability  of  the  many  beamlets  which  rapidly  merge  into  a well  defined  beam.* 

A careful  observation  of  the  beam  structure  reveals  that  the  waist  of  the  beam 
is  fairly  constant  near  the  downstream  face  of  the  accelerator  and  for  a few 


* See  on  this  subject  Refs.  26,  27,  65,  68,  71. 
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diameters  away,  and  that  the  ion  stream- lines  bend  slightly  outward  to  reach  a 
fairly  constant  expansion.  Therefore,  the  virtual  origin  of  the  conical  envelope 
is  shifted  forward  while  the  actual  stream-lines  are  slightly  curved  rather  than 
ideally  straight.  Closer  in,  therefore,  the  model  would  overestimate  the  maximum 
transverse  velocity  component  accessible  to  the  ions  and  thus,  according  to  the 
synthetic  temperature  scheme,  the  speed  ratio  simulated  would  be  underestimated 
as  indicated  by  the  measurements.  The  model  is  recovered  at  about  z/r0  ~ 30. 
Conversely  because  of  the  curvature  of  the  ion  stream- lines  the  transverse 
velocity  components  would  be  underestimated  beyond  this  station,  thus  leading 
to  an  over-estimation  of  the  speed  ratio  which  is  also  observed.  At  large 
axial  distances,  however  (z/rG  S'  90),  the  speed  ratios  measured  using  the  finite 
probe  effect  seem  to  level  off.  This  also  is  believed  to  be  a real  flow-field 
effect  rather  than  a lack  of  refinement  of  the  method.  It  is  typically  expected 
for  speed  ratios  larger  than  roughly  100,  for  probes  with  an  aspect  ratio  (i/d) 
greater  than  about  100,  as  indicated  on  Figs.  1-7  and  1-8  where  the  nondimensional 
probe  length  rg  is  near  unity.  Such  an  indication  was  also  suggested  by  an 
experimental  investigation  reported  in  the  literature  (Ref.  78)  which  showed  that 
the  ion  Mach  number  measured  by  scanning  cup- analyzers  levelled  off  to  a ’saturated’ 
value  for  large  axial  distances  downstream.  The  ion  source  was  similar  in  concep- 
tion to  the  one  used  for  this  experiment.  The  weight  of  evidence  at  this  stage  is 
to  ascribe  the  measured  results  to  a non-perfect  flow  model  rather  than  to  a defect 
in  the  probe  theory.  To  shed  more  light  on  this  we  will  present  a more  complete 
correlation  of  the  results. 


IV. 3.2  Reliability  of  the  Method  to  Measure  the  Speed  Ratio 
(a)  Ion  Source  Effects 


In  order  to  account  for  any  spurious  effect  originating  in  the  ion 
source  itself  which  could  perturb  the  probe  signal,  we  have  grouped  a set  of 
data  points  obtained  for  a value  of  z/r0  of  40  for  different  beam  energies. 
The  ion  speed  ratios  which  are  obtained  from  the  self-consistent  solution  at 
a given  ion  beam  energy  are  normalized  for  convenience  with  respect  to  the 
speed  ratio  measured  at  a reference  beam  energy  of  700  volts.  Such  a corre- 
lation of  the  results  is  presented  in  Fig.  IV-l4.  The  choice  of  S/rQ  = 40 
provides  an  overlapping  of  the  results  obtained  with  the  two  grid  sets  which 
were  used  to  cover  a wider  range  of  z/rQ  by  merely  doubling  the  beam  diameter: 


2rQ  =2.54  cm; 

7 ~ 10’ 

2r  = 5.08  cm; 
0 ’ 

7 ~ l4‘ 

A slight  variation  of  the  beam  divergence  angle  7 has  not  been  thought  to 
introduce  any  significant  variation  in  the  measured  values  of  Si  for  the  two 
beam  sizes.  It  is,  however,  evident  from  the  trend  of  Si  that  the  ion  source 
does  introduce  a spurious  beam  energy  dependence  for  energies  lower  than  about 
300  eV.  Such  an  artificial  dependence  of  Si  on  beam  energy  must  therefore  be 
kept  in  mind  when  attempting  to  properly  scale  the  ion  temperature  in  the  flow- 
field,  especially  for  simulation  purposes.  For  beam  energies  larger  than  300 
eV,  the  results  are  nevertheless  consistent  with  one  another  within  about  15$ 
up  to  the  maximum  value  of  beam  energy  we  have  used  (1200  eV)  . 
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When  one  considers  the  relatively  complex  mechanisms  responsible  for 
the  emission  of  the  ions  into  a well  defined  beam  (see  Section  III. 1.2),  it  is 
not  surprising  to  obtain  a certain  modification  of  the  beam  structure  for  the 
lower  part  of  the  accelerator  voltage  range.  In  fact,  it  is  well  known  from 
experience  with  ion  engines  that  the  beam  formation  becomes  imperfect  in  the 
lower  extraction  voltages.  This  behaviour  can  be  demonstrated  further  when  the 
low-voltage  accelerator  configuration  is  used.  In  this  mode  the  optics  of  the 
beam  forming  apertures  are  less  constrained  because  there  is  only  a single 
extraction  grid  (Fig.  II-5) . The  ion  speed  ratios  measured  under  these  conditions 
are  presented  in  Fig.  IV-15  where  the  probes  were  again  positioned  at  a z/r0  value 
of  40.  A stronger  net  beam  voltage  dependence  is  shown  in  this  configuration. 

(b)  Probe  Bias  Effects 


We  can  further  estimate  the  reliability  of  the  method  to  measure  Si 
by  noting  the  effect  of  a substantial  variation  of  the  probe  bias  on  the  values 
of  Sj_  obtained  from  the  probe  response.  Three  sets  of  data  points  are  presented 
on  Fig.  IV-l6  which  show  that  the  probe  bias  does  not  affect  the  speed  ratios 
generated  by  the  method  significantly.  No  systematic  influence  is  observed  for 
a probe  bias  range  which  covers  in  excess  of  one  order  of  magnitude  change  in 


0.03  < <t>  <0.5 
P 

We  have  added  three  experimental  points  obtained  from  Ref.  57  for  comparison, 
in  which  Fournier  and  Pigache  report  some  attempts  to  verify  Sanmartin's  finite 
probe  theory.  These  quantitative  attempts,  though  limited  in  number  were  made 
to  obtain  an  estimate  of  the  ion  speed  ratio  according  to  Sanmiartin' s analytical 
description.  The  authors  did  not  use  a self-consistent  scheme  but  actually 
estimated  the  value  of  6 which  they  introduced  in  the  analytical  expressions.  The 
significant  dependence  on  <Bp  displayed  in  Fig.  IV-16  by  their  data  points  could 
possibly  be  tied  to  an  incorrect  choice  for  "5  (see  Section  IV .4 .2). 

(c)  Probe  Aspect  Ratio  Dependence 

To  complete  our  test  of  the  degree  of  refinement  obtained  by  the  method 
proposed  to  determine  the  ion  speed  ratio  using  the  end  effect,  we  must  verify  the 
correlation  between  the  values  of  Sj_  determined  by  the  use  of  significantly 
different  probe  aspect  ratios.  This  information  was  already  available  in  Fig. 

IV- 13.  The  variations  between  the  values  of  Sj_  measured  by  six  different  probe 
i/d  ratios  are  repeated  in  Fig.  IV-17.  These  values  of  are  obtained  at  a 
fixed  z/rQ  which  is  equivalent,  according  to  our  geometrical  model,  to  keeping 
Sj_  constant.  These  typical  results  indicate  that  the  values  of  Si  determined  by 
direct  measurement  of  the  probe  angular  response  are  consistent  within  about  20 $> 
for  a probe  aspect  ratio  greater  than  about  70.  The  anomalous  behaviour  for  the 
very  low  aspect  ratio  probe  is  again  quite  different. 

Using  the  same  group  of  data  points  we  can  also  illustrate  the  fact 
that  the  method  leads  to  a reasonable  agreement  between  the  values  of  Si  detected 
by  different  probe  sizes.  Figure  IV-18  which  compares  measurements  at  two  dif- 
ferent stations  (or  nominal  speed  ratio  values)  suggests  that  the  absolute  probe 
diameter  does  not  introduce  any  systematic  dependence  on  the  measurement  of  Si, 
at  least  for  the  range  studied  here.  In  particular,  the  group  of  points  at  4,  8 
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and  10  mils  which  correspond  to  probes  with  roughly  the  same  probe  aspect  ratio 
(i/d  ~ 170)  was  purposely  selected  to  corroborate  this  fact.  The  data  points 
shown  at  20  mils  correspond  to  the  anomalous  R ~ 22  case,  therefore  no  conclusion 
about  the  influence  of  the  probe  diameter  on  the  integrity  of  the  method  can  be 
drawn  from  this  particular  probe.  Since  the  importance  of  the  end  effect  decreases 
with  increasing  probe  diameter,  at  a given  set  of  stream  conditions,  one  can  expect 
the  sensitivity  of  the  finite  probe  to  Sj_  to  decrease  accordingly  (see  Fig.  1-8). 
However,  we  can  estimate  that  a 20  mil-diameter  probe  with  a value  of  i/d  of  170 
would  not  have  yielded  a measurement  of  Sj_  significantly  different  from  the  ones 
obtained  with  the  probes  used  here. 


IV. 3. 3 Reconstruction  of  the  Probe  Response 

In  order  to  implement  Sanmartin's  theoretical  analysis  of  the  end 
effect  and  to  demonstrate  the  overall  ability  of  the  self-consistent  method  which 
we  developed  to  obtain  the  ion  speed  ratio  from  the  angular  response  of  cylindrical 
probes,  we  have  generated  a simple  computer  code  which  reconstructs  the  entire 
probe  response.  It  requires_the  values  of  the  beam  and  of  the  probe  parameters, 
together  with  the  values  of  5 and  of  Sq  obtained  by  the  self-consistent  procedure. 
Thus,  by  using  the  following  features  measured  from  the  experimental  traces, 

j(e  = o),  Jje  = 0),  J(9  = 77/2),^ 

the  entire  current  angular  response,  i.e.  over  the  entire  angular  range, 


0°  < 0 < 90* 


can  be  reproduced  as  shown  in  Fig.  IV-19.  The  computer  plots  are  scaled  according 
to  the  real  current  scale  such  that  the  normalizing  value  of  = 0)  corresponds 

to  unity.  Two  curves  are  plotted  for  each  comparison.  First  the  peak  structure 
is  traced  for 


0 < 0 < 30* 


and  then  the  'infinite'  probe  response  is  superposed  on  the  graph  for 


0°  < 0 < 90° 


in  order  to  check  the  extrapolated  region  as  shown  on  the  lab  trace.  A scale 
amplification  factor  of  4 is  also  used  for  the  angle  axis  to  compare  the  theory 
with  the  experimental  traces  which  were  also  taken  with  an  identical  angular 
sensitivity  (x  4)  in  order  to  measure  the  behaviour  at  small  angles  more 
accurately. 

A set  of  12  cases  were  selected  for  comparison.  This  set  covers  the 
complete  range  of  the  parameters  studied  and  is  listed  in  Fig.  IV-20  with  the 
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corresponding  parameters.  The  direct  comparisons  are  shown  on  Figs.  IV-2l(a)  to 
(e)  where,  the  theoretical  curves  are  matched  to  the  measured  traces.  Some  extra 
features  are  added  such  as  a thick  vertical  bar  on  the  angle-axis  locating  the 
predicted  value  for  the  limit  to  the  range  of  validity  of  the  theory  and  a thin 
horizontal  line  locating  the  half-widths  according  to  the  nomenclature  introduced 
in  Fig.  IV- 1.  The  open  symbols  superposed  on  the  theoretical  curves  are  discrete 
measurements  of  the  current  response.  In  each  case  the  actual  experimental  curve 
is  shown  for  one  half  of  the  complete  angular  range  from  -90*  to  +90*.  To  obtain 
a clear  exposition  the  other  half  of  the  angular  range  shows  the  theoretical  result 
as  a-  solid  curve,  with  superposed  discrete  experimental  measurements.  Each  case 
includes  the  normal  and  the  X4  blow-up,  for  both  theory  and  measurement.  After 
inspection  we  can  conclude  that,  in  general,  the  agreement  is  excellent. 


We  can  see  that  the  theory  indicates  that  the  peak  current  approaches 
zero  for  large  angles  and  thus  it  fails  to  describe  the  peak  structure  correctly 
for  values  of  pitch  angle  which  are  beyond  about  two  half-widths  from  the  maximum 
(0  =0).  If  we  take  the  more  natural  definition  for  the  limit  to  the  validity  of 
the  theory  as  the  angle  at  which  the  peak  prediction  equals  the '.infinite'  probe 
value,  we  make  the  following  observation  from  the  12  typical  cases  depicted  in 
Fig.  IV-21  - namely,  that  the  theoretical  limit  obtained  by  Sanmartin, 


as  indicated  by  the  thick  vertical  bar,  systematically  overestimates  the  inter- 
secting points  of  the  two  calculated  curves.  In  fact  we  have  obtained  the 
corresponding  intersection  experimentally  by  smoothly  extrapolating  both  the 
decrease  in  peak  current  and  the  'infinite'  probe  current  on  each  lab  trace.  The 
comparisons  for  the  location  of  this  intersection  are  shown  in  Fig.  TV-22  where 
the  full  and  the  open  symbols  are  the  intercepts  obtained  for  the  measured  and  the 
calculated  casesj  respectively.  .Since  in  the  expression  for  0]_im  the  leading  term, 
proportional  to  6,  dominates  for  values  of  > 20,  this  curve  appears  to  favour 
values  of  5 smaller  than  those  calculated  using  Sanmartin' s approximation  (Eq. 
1.39)  by  about  40$.  A similar  conclusion  results  from  a direct  comparison  of  the 
values  of  6 obtained  from  the . self-consistent  calculational  process  with  the  values 
from  Sanmartin' s approximation. 


A relation  between  the  nondimensional  probe  bias  Gt,  and  & can  be  obtained 
using  Eq.  I. 39. 


$ ~ exp 
P 1 


25_1  - in  5 + 0.307  - 2in 


(IV. 7) 


This  relation  is  plotted  in  Fig.  IV-23,  with  6*  as  a parameter,  together  with  a 
number  of  experimental  values  of  <£p_,  &Self  and.  €*.  These  data  points  are  referring 
to  the  self-consistent  values  for  5 obtained  using  our  regular  procedure.  A 
direct  comparison  using  the  full  circles  for  instance,  shows  that  for  any  of  the 
given  <£p  reported  here  for  an  e*  ~»01  the  value  of  the  self-consistent  6 which  is 
obtained  is  smaller  than  the  corresponding  approximated  value  derived  from  the 
curve  by  about  30$.  This  disagreement  is  believed  to  be  due  to  a readjustment 
necessitated  by  the  fact  that  the  real  flow  field  deviates  from  the  assumed 
conical  cum  Maxwellian  models  and  thus  requires  that  the  theoretical  value  be 
slightly  modified. 


Non-Maxwellian  Effects 


An  effect  which  is  not  accounted  for  by  the  theory  is  materialized  in 
the  occurrence  of  side-peaks  on  the  probe  response.  These  peaks  have  a much 
smaller  amplitude  than  the  main  peak.  They  actually  could  only  be  detected  with 
the  long  slender  cylindrical  probe  (Fig.  IV-24)  or  with  the  small  collimated 
Faraday  cups  (Fig.  A-5).  Such  secondary  peaks  were  also  observed  by  Fournier 
and  Pigache  (Ref.  57)  on  the  current  angular  response  of  their  small  cups.  In 
order  to  reconcile  large  discrepancies  between  the  speed  ratios  they  estimated 
by  a method  briefly  outlined  in  our  Appendix  D,  and  some  rough  estimates  they 
obtained  using  Sanmartin's  analysis,  they  first  alluded  to  the  possibility  of 
non-Maxwellian  effects.  Indeed,  Sanmartin's  analysis  was  based  on  the  assumption 
of  an  ideal  thermal  ion  population  (drifting  Maxwellian)  far  ahead  of  the  probe. 
The  concept  of  thermalization  introduced  in  our  stream  model,  however,  rests  on 
an  artificial  basis.  The  synthetic  ion  population  obtained  from  the  model  we  have 
adopted  is  likely  to  be  described  with  a distribution  which  is  different  from  a 
Maxwellian  distribution.  The  effect  introduced  is  relatively  small.  Thus  we  are 
led  to  believe  that,  since  the  results  are  not  basically  in  error,  the  deviation 
from  ideal  Maxwellian  does  not  lead  to  serious  conpli cations  and  furthermore  that 
the  use  of  the  self-consistent  method  appears  to  provide  some  of  the  required 
adjustments  automatically. 


The  fact  that  two  very  much  different  probe  structures  wouud  detect 
the  secondary  peaks  in  a similar  fashion  is  not  so  surprising.  Indeed,  both 
ion  collection  schemes  are  sensitive  to  the  transverse  component  of  the  ion 
population.  Therefore,  these  observations  must  reflect  some  anomaly  occurring 
in  the  flow  itself  since  they  are  not  restricted  to  a particular  probe  type. 


Although  only  limited  quantitative  arguments  can  be  invoked,  we  can 
attenpt  to  elucidate  the  question  of  the  origin  of  the  secondary  peaks  by  locating 
at  what  angles  the  measured  side-peak  features  are  occurring  on  the  traces.  These 
angles  are  indicated  in  Figs.  IV-24  and  A-5*  We  notice  that  they  are  located  at 
angular  positions  which  are  beyond  the  angular  width  of  the  source  exit  diameter  D0 
as  sighted  at  the  given  probe  location  (typically  2°  to  5°)*  These  positions  of 
the  side-peaks  correspond  to  angles  which  are  close  to  the  spreading  angle  y of 
the  source  in  all  the  cases  reported  here.  We  have  also  reported  in  the  table  below 
some  results  obtained  from  Ref.  57  which  also  show  secondary  peaks  although  under 
quite  different  experimental  conditions.  These  results  are  summarized  in  the 
following  table: 


Long  Cylindrical  Probe  (Fig 


acondary  Peak  Location 


Collimated  Faraday-Cup  Probe  (Fi 


31 

42 

182  (Ref .57) 

~ 8.2°  (vB  = 500V) 
~ 9*2°  (vB  = 500V) 
~ 3.5-4°  (vB  ~ 20V) 

7 - 10°  (D0  =2.5  cm) 
7 ~ 10°  (D0  = 2.5  cm) 
7 ~ 10°  (D0  = 10  cm) 

We  notice  finally  that  the  side-peaks  are  vanishingly  small  for  the  larger  beam 
energies  as  shown  in  Fig.  IV-24.  This  effect  likely  originates  at  the  accelerator 
grids.  Although  no  photographic  records  have  supported  this  possibility,  it  is 
conceivable  that  a small  sub-beam,  fringing  symmetrically  outward  within  the  main 
beam  is  produced  because  of  a fortuitous  inhomogeneous  ion  density  at  the  grid  and 
that  the  effect  is  compounded  at  lower  beam  energies  (see  Section  IV. 3 .4a).  In 
final  analysis,  however,  the  observation  of  these  small  secondary  features  is  an 
indication  of  the  capability  to  resolve  find  structure  with  the  probe  effect. 


IV. 3-5  Estimate  of  Errors 


Our  confidence  in  the  ion  speed  ratios  measured  using  the  ion  temperature 
sensitivity  of  the  finite  cylindrical  probe  is  supported  by  the  fact  that  the 
results  were  found  to  be  consistent  within  our  conservative  estimate  of  the  margin 
of  error  even  though  some  departure  from  our  simplified  model  was  found  to  be 
significant  for  the  larger  values  of  Si . 

In  order  to  get  this  estimate  of  the  margin  of  error,  we  can  refer  to 
the  schematic  representation  of  the  method  we  have  used  to  obtain  both  the  charac- 
teristic field  length  6,  and  the  ion  speed  ratio  Si,  according  to  a self-consistent 
scheme  as  shown  in  Figs.  IV-8  and  TV-9.  The  arguments  leading  to  this  estimate  are 
thus  intimately  related  to  this  direct  formulation  of  the  solution.  One  can  see 
that  the  iterative  procedure  yielding  the  self-consistent  values  for  6 and  Si  is 
hinged  firstly,  on  the  extrapolation  of  the  'infinite’  probe  current  value  at 
9 = 0,  secondly  on  the  evaluation  of  the  slow-ion  offset  which  by  the  same  token 
sets  the  true  half-height  and  thereby  locates  the  corresponding  peak  half-width, 
and  thirdly  on  the  agreement  between  the  peak- shape  functions  Xi  and  X2.  After 
a closer  inspection  of  Fig.  IV-9  one  notes  that  the  determination  of  the  peak- shape 
function  X2  obtained  by  the  relation  involving  the  measured  half-width  01  and  given 
by  2 


x2 


(IV. 8) 


is  not  strongly  dependent  on  5 for  the  range  of  interest,  6 > 0.1,  whereas_Xi 
obtained  via  the  peak  maximum  (0  = 0)  using  Eq.  1,26  is  more  sensitive  to  8. 

This  is  reflected  in  the  expression  generated  in  order  to  represent  the  relative 
error  in  both  Xi  and  X2  using  the  analytical  representations  provided  by  1*26, 
1,29  and  IVk8.  We  get 


AXi  „ AR  1 / \ 

Xl  - I + 2 \ -IT) 


+ 


(iv. 9) 


2 


(iv.10) 
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where  the  symbol  0 indicates  an  error  of  second  order.  Now  since  the  iterative 
sequence  is  generated  by  comparing  Xi  and  X2  in  turn,  it  is  thus. reasonable  to 
use  Eas.  IV. 9 and  IV. 10  to  yield  an  effective  relative  error  on  5 by  putting 


The  relative  error.on  R is,  in  fact,  not  likely  to  introduce  a significant 
contribution  in  (A6/5)  since  it  would  simply  account  for  a variation  acting  in 
the  'same  direction'  in  both  Xi  and  X2  whereas  the  residual  error  *0  could  be 
contributing  in  'opposite  directions'  and  thus  the  worst  case  is  chosen. 


Similarly  an  expression  for  the  relative  error  in  S^  can  be  arrived  at. 
First  we  rearrange  Eq.  1.28,  i.e.,  we  get 


S.  = — 
1 ^ 


2R 


L (m  + 0.6Tg2)ilncr2  J 


Moreover  for  the  experimental  range  considered  here  we  can  estimate,  using  the 
analysis  provided  in  Section  1.1.2,  that 


X cir 


(IV.12) 


Therefore  the  relative  error  in  X can  be  assessed  directly  since  it  will  contribute 
to  the  error  in  Sj_  via  the  final  calculation  of  S^ , as  shown  in  the  flow-chart 
depicted  in  Fig.  IV-8,  i.e., 


x - 2 n u 

which  means,  using  the  expression  obtained  above  for  Si(t]),  for  the  range  of 
interest  Tg  < 1,  we  have 


i .e 


5 


Si"2V  T| 


/V/ 

+ 0 


~ to  , ar 

Si  “ X R 


+ 0 


Once  the  value  of  5se^f  has  been  established  by  iteration,  it  follows  that  a 
conservative  assessment  of  the  relative  error  in  S^  can  be  attempted  using  the 
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expression, 


Finally,  we  can  obtain  an  estimate  of  the  relative  error  in  each  of  the 
parameters  contributing  to  A6/S  and.  AS 

(1)  The  residual  error  contained  in  0 is-  not  expected  to  exceed  4$,  according  to 
the  analytical  treatment,  and  is  averaged  at  about  2$. 

(2)  Using  systematic  measurements  of  the  various  cylindrical  probes  the  average 
error  on  the  probe  aspect  ratio  is  evaluated  using 


AR  _ A0  Ad  ~ * 

(3)  The  relative  error  in  the  determination  of  the  probe  bias  can  be  calculated 
from  the  relation  we  obtained  in  Chapter  III,  i.e., 

V " ve 

<D E l 

P " VA  + - V, 


Since  we  use  V ~ 0,  and  moreover  we  have 
P 


|VD  - Vf  I « v+ 


then 


A$>  Av„ 

~ f 


AV. 


J,  - Vf+  v^<6ii 


upon  using  the  standard  meter  reading  accuracies. 

(4)  The  current  peak  angular  accuracy  was  evaluated  at  less  than  about  4 % which 
corresponds  to  an  average  absolute  error  of  less  than  0.3°;  thus 


(5)  The  relative  error  in  the  non dimensional  current  peak  height  j/j^  (0  = 0)  is 
evaluated  using  the  relation  which  was  obtained  for  the  determination  of  the 
background  slow-ion  contribution 


(J(o)  - J) 
(Jjo)-  J) 


( IV. 14) 
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, ~ ™ — 


where  the  subscript  (N)  refers  to  the  ideal  conditions 
background  offset  is  present.  A conservative  estimate 
treatment  of  the  current  traces  and  typical  conditions 
yields 


A(  ^ (0) 

00 

f <°> 

00 


< 


prevailing  when  no 
using  our  systematic 
was  obtained  and 


This  relatively  large  value  is  due  to  the  fact  that  this  quantity  is  obtained 
by  direct  extrapolation  on  the  trace,  and  that  the  slow-ion  offset  must  be 
estimated  using  this  extrapolated  value  of  J^O) . The  situation  is  complicated 
in  cases  where  J^O)  is  a relatively  small  value,  comparable  say  to  J,  such 
that  the  relative  error 


A[JjO)  - J] 

j (o)  - IT ' 

00 x 7 

can  be  large.  In  such  cases  multiple  readings  are  obtained  to  improve  the 
situation. 


In  conclusion  the  net  relative  error  on  the  estimate  of  Sq  using  our 
method  can  be  set  at  less  than  about  32$.  Consequently  the  method  is  expected 
to  yield  ion  temperatures  with  an  accuracy  in  the  neighbourhood  of  60$.  This  is 
not  such  a prohibitively  large  error  when  one  considers  that  the  end  effect  could 
very  well  be  the  only,  if  not  the  simplest  probe  characteristic  response  from 
which  the  ion  temperature  can  be  directly  extracted.  In  fact  the  weight  of  the 
evidence  presented  in  the  various  discussions  carried  out  in  this  last  section 
shows  that  the  measured  values  of  Si  were,  in  general,  consistent  with  one  another 
with  an  accuracy  better  than  20$.  This  indication  is  supported  by  the  fact  that 
the  conservative  error  estimates  considered  here  were  'worst  case'  estimates  in 
the  sense  that  no  statistics  were  introduced  in  the  discussions.  This  is  especially 
true  for  the  measurement  of  the  peak  features  as  described  in  Section  IV. 1.1.  As 
indicated  for  the  various  typical  cases  presented  in  Figs.  IV-2l(a)  to^(e) , we 
have  obtained  the  extrapolated  value  of  9 = 0)  and  the  half-width  01  from  both 
the  'standard'  trace  and  from  the  'blown-up'  trace  in  each  case,  and  from  repeated 
traces  in  a number  of  cases,  in  order  to  improve  the  sampling  of  the  measurements. 


It  can  also  be  expected  that  the  level  of  confidence  in  the  measurement 
of  5,  which  is  also  allowed  by  our  method,  could  be  somewhat  higher  than  the 
estimate  suggests,  such  that  the  determination  of  this  important  field  parameter 
could  be  obtained  with  an  accuracy  comparable  with  the  one  suggested  by  the 
experimental  results,  namely  20$.  Finally  it  was  found  that  the  operation  of 
the  ion  source  in  a testing  facility  with  as  low  as  possible  a 'background'  ion 
content  would  be  highly  desirable  in  order  to  conduct  a diagnostic  of  the  kind 
reported  here.  In  final  analysis  we  have  seen  that  the  bulk  of  the  uncertainty 
in  Si  is  introduced  when  an  offset  current  due  to  charge-transfer  mechanisms  is 
superposed  on  the  sampled  fast-ion  current  to  the  probes.  Nevertheless  one  could 
obtain  a reasonably  simple  scheme  to  account  for  this  offset,  and  to  estimate  it 
off  the  experimental  traces,  without  jeopardizing  the  overall  reliability  of  a 
diagnostic  method,  such  as  the  one  we  have  adopted  in  this  investigation  to 
measure  S. . 
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CONCLUSIONS 


An  experimental  facility  providing  a flexible  range  of  high-speed, 
low-density  plasma  conditions  has  been  designed  and  constructed  for  a number 
of  applications  ranging  from  the  simulation  of  flight  in  the  ionospheric 
environment  to  the  proving  of  diagnostic  techniques  for  various  laboratory 
and  industrial  plasmas.  The  operation  of  the  plasma  source  is  reliable  and 
requires  only  limited  ancillary  instrumentation. 

The  end  effect  theory  proposed  by  J.  R.  Sanmartin  has  been  validated 
and  the  various  parametric  dependences  have  been  determined  experimentally.  A 
wide  range  of  flow  conditions  and  of  probe  dimensions  were  considered.  In 
particular,  this  experimental  verification  was  obtained  over  a considerably 
wider  range  of  probe  aspect  ratios  ($/d)  and  of  ion  speed  ratios  than  was 
previously  investigated.  Moreover  the  effect  was  studied  for  significantly 
warmer  ions  (T^/Te  ~ l)  whereas  all  the  cases  reported  until  now  had  been 
obtained  for  basically  cold  ions  (Tq/Te  ~ 0) . Sanmartin' s theory  was  described 
in  this  report  in  terms  of  a new  set  of  parameters  more  readily  accessible  from 
direct  physical  measurements,  and  which  demonstrates  explicitly  the  relatively 
weak,  sensitivity  of  the  effect  to  the  electron  temperature.  The  limits  to  the 
range  of  validity  of  the  theory  have  been  outlined  using  this  new  set  of  para- 
meters . 


A method  to  measure  the  ion  temperature,  directly  from  the  angular 
response  of  a cylindrical  probe,  in  a collisionless  plasma  stream,  has  been 
demonstrated.  The  height  and  the  half-width  of  the  ion  peak,  which  are  both 
sensitive  to  the  ion  speed  ratio,  are  used  simultaneously  to  ascertain  that 
value  in  a self-consistent  way.  The  experimental  results  were  checked  against 
a simple  geometrical  model.  This  has  led  us  to  conclude  that  the  flow-field 
model  proposed  was  not  correct  as  the  ion  streamlines  were  found  to  bend  slightly 
downstream  of  the  source,  while  the  model  assumed  a single  conical  expansion. 

This  evidence  was  further  corroborated  after  inspection  of  a series  of  radial 
beam  profiles.  However,  a cross- correlation  of  the  results  indicated  that  the 
measurements  of  the  ion  speed  ratio,  obtained  by  various  probes  and  under  many 
flow  conditions,  were  consistent  with  each  other  within  an  acceptable  margin 
of  about  20%.  In  fact  the  usefulness  of  the  probe  would  be  expected  to  be 
substantially  improved  under  the  ionospheric  conditions  for  instance,  since  the 
measurements  would  not  suffer  from  the  inherent  contribution  of  spurious  slow- 
ions  on  the  current  collection,  as  it  is  the  case  with  extended  testing  tunnels. 
Moreover  the  successful  theoretical  reproduction  of  the  entire  structure  of  the 
probe  angular  response,  which  has  been  derived  for  a number  of  typical  cases, 
epitomized  the  overall  reliability  of  the  method.  Finally  we  have  observed 
that  the  long  slender  probes  were  able  to  detect  interesting  'fine-structure' 
aspects  of  the  drifting  ion  population;  indeed  the  small  secondary  ion  peaks 
detected  on  the  probe  angular  response  suggested  that  non- Maxwellian  effects 
were  introduced  in  the  beam  prior  to  the  current  collection. 

Some  areas  of  applications  of  the  end  effect  can  therefore  be  outlined. 
As  a typical  case,  the  method  can  be  extended  to  multiple  gas  species  in  order 
to  determine  the  ion  temperature  of  the  heavy  constituents  of  the  ionosphere  for 
altitudes  higher  than  about  300  km  (H+,  He+,  N+,  0+...).  This  could  have  a 
significant  advantage  since  the  strong  attituae- sensitivity,  easily  identified 
by  the  sharp  ion  peaks,  could  be  proved  useful  as  a sensing-and-control  feature 
in  order  to  obtain  the  fine-pointing  accuracy  required  by  sophisticated  instrument 
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packages.  This  sensitivity  to  attitude  and  to  random  motion  within  the  plasma 
stream  can  also  be  useful  to  the  development  of  sources  such  as  ion  thrusters. 

In  fact,  these  devices  which  had  been  developed  exclusively  for  propulsion  missions 
in  the  mid-sixties  (Refs.  26,27)  are  now  recognized  as  outstanding  candidates  for 
a number  of  industrial  applications  requiring  an  advanced  quality  control,  such 
as  etching  and  thin- film  deposition  techniques  (Ref.  87),  or  even  as  large-scale 
plasma  injectors  for  fusion  machines  (Refs.  85,86). 


It  is  suggested  here  that  a number  of  interesting  aspects  of  the 
observation  of  non-Maxwellian  effects  on  'standard'  probe  characteristics  could  be 
elucidated  by  introducing  a non-Maxwellian  treatment  in  Sanmartin's  analysis.  This 
could  eventually  lead  to  the  correct  description  of  other  probe  'distortions'  under 
non-ideal  mass-motion  conditions.  Finally,  as  an  interesting  and  stimulating 
extension  of  this  report,  the  experimental  validation  of  Sanmartin's  analysis 
reported  in  Ref.  23  should  be  carried  out  for  the  case  of  a quiescent  collisionless 
plasma.  Indeed,  the  analytical  description  of  the  ion  temperature  sensitivity  of 
the  end  effect  allows  the  conceptual  equivalence  between  the  steady-flow  problem 
and  the  one  dimensional  time- dependent  problem  wherein  an  'infinitely'  long 
negatively  biased  cylinder  is  immersed  instantaneously  in  an  unperturbed  plasma 
at  rest  and  the  now  familiar  ion  current  overshoots  are  observed. 
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APPENDIX  A - MASS-FLOW  MONITORING 


A simple  and  reliable  system  was  constructed  to  control  and  monitor 
the  very  small  argon  mass-flows  required  for  this  experiment  (Z  0.1  ni/sec) . 
The  system  is  composed  of  a known  volume  V with  a source  exit  pressure 

Pg  « Py  (ideql  pump) 

and  the  feed- line  volume  is  negligible.  The  system  temperature  is  maintained 

constant  at  T , thus 
o 

P = P 

° (at  t = 0) 

T = T 

o 

(l)  The  mass  flow  at  steady-state  is  given  by  the  gas-law  relationship 


W = 


W + _ W f V ~ 

RT  RT  V “ At  ) 


(2)  With  T = 300°K  and  argon  gas  from  a volume  V = 5.175,0,  we  have 


W ~ 8.5  x 10~2  P (ipg/sec) 

for  £ (Nt/n^/sec)  or 

W ~ 0.45  P (equiv.  anpere)* 

with  P (torr/min) 

(3)  Therefore,  the  mass-flow  is  known  by  monitoring  the  pressure  drop  of  the 
reference  volume  with  time,  provided  the  total  drop  is  small  with  respect  to 
the  initial  pressure  (see  Fig.  A-l) 

AP  « P 

o 

Typically,  with 

PQ  ~ 750  torr 

At  ~ 5 hours 

^ (l  torr/20  min) 

then,  we  verify  (3)  since 


p 

o 


Zf, 


* 


The  equivalent  airpere  is  a quantity  which  assumes  that  6.25 
atoms  represent  as  many  single-charge  units. 


x 


10  1 8 neutral 
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For  argon  we  have  the  relations 

1 eq.  amp  = 0.42  mg/sec 
1 eq.  amp  = 0.187  torr-i/sec  @ STP 
1 eq.  amp  = 0.178  |ii/sec  @ STP 
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OUTLINE  OF  THE  METHOD  DEPICTING  THE 
APPROXIMATE  SOLUTION  of  (Y,  - Y,  ) YIELDING 


length  ratio) 


a)  ESTIMATED  values  or  £ ,n,,Te  using  the  _ 
ELECTRON  TEMPERATURE  sensitivity  from  6 


APPENDIX  B - PROBE  SPURIOUS  EFFECTS 


Both  magnetic  and  R.F.  fields  can  he  neglected  here;  the  spurious 
effects  on  the  collection  of  positive  ions  can  be  described  briefly  under  three 
headings : 


- ion  current  over- estimation 

- ion  current  under-estimation 

- contamination 

(i)  Ion  Current  Over-Estimation 

Over-estimation  may  originate  from  secondary-electron  emission  from 
the  probe  due  to  thermionic  effects,  photo-emission,  ion  and  metastable  atom 
collisions  or  ion  and  atom  bombardment  (2nd  Townsend  yield)  at  the  metallic 
surface  (Ref.  79)-  The  first  three  sources  can  be  neglected,  from  rate-process 
probabilities  and  surface  energy  considerations  (Ref.  13),  though  ion  and 
metastable  atom  bombardment  is  occurring  in  this  regime.  Secondary-electron 
yield  from  direct  bombardment  would  constitute  the  major  contribution.  Reported 
figures  indicate  a maximum  yield  of  about  0.1  electron/ion  at  incident  energies 
up  to  1.5  KeV,  for  singly  ionized  argon  on  atomically  clean  W.  Gas  covered 
('dirty')  conditions,  which  occur  when  foreign  atoms  are  adsorbed  on  the  surface, 
would  give  a figure  generally  much  lower  (Ref.  63)  and  an  estimate  of  0.05  is 
representative  of  the  maximum  yield  in  this  regime.  Similarly  met as table- atom 
bombardment  would  not  imply  any  significant  contribution.  From  these  conservative 
estimates  therefore,  no  correction  was  made  to  the  collected  currents.  Finally, 
more  insidious  effects  such  as  reflections,  oscillations  and  'trapped'  orbits 
have  been  neglected. 

(ii)  Ion  Current  Under-Estimation 


Since  the  probe  is  made  strongly  negative  for  the  present  study,  this 
second  category  is  more  likely  to  be  less  important  than  the  sources  indicated 
under  (i).  Ion  reflections  or  ejections  and  sputtering  of  heavy  particles 
(generally  neutrals  from  the  metal)  although  significant  in  some  instances 
(see  Appendix  C),  is  of  little  concern  to  us  here.  Finally,  from  conservative 
estimates,  doubly- charged  ions  would  not  contribute  to  more  than  about  4$  of 
the  net  currents. 

(iii)  Contamination 

This  subject  has  been  treated  in  Appendix  C,  yet  it  is  worth  noting 
that  contamination  of  the  probe  can  account  for  some  error  since  it  can  imply 
that  the  true  collecting  surface  is  not  known  precisely.  Conducting  or  insulating 
layers  on  the  collector,  or  on  the  ceramic  stem,  could  originate  from  a few  sources 
backstreaming  oil,  evaporation  of  tungsten  from  the  source  and  the  neutralizer 
filaments,  and  possibly  from  backspottered  material  (copper  and  aluminum)  from 
within  the  test  section;  the  knowledge  of  such  possibilities  has  led  us  to  regular 
and  careful  inspection  of  the  probes.  A black  layer  of  baked  oil  was  occasionally 
observed  on  the  insulator  parts  and  found  to  be  slightly  conductive.  The  probe 
was  then  flame  cleaned  and  fresh  ceramic  cement  was  applied.  Both  hot  tungsten 
emitters  required  on  the  accelerator  were  screened  to  prevent  direct  line- of- sight 
with  the  probe.  No  significant  metallic  coating  was  observed  for  extensive  running 
periods . 
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APPENDIX  C - ION-BOMBARDMENT  PROBE  CLEANING 


The  subject  of  probe  cleaning  for  diagnostic  purposes  is  a difficult 
subject  on  which  to  obtain  full  agreement  from  all  probe  users.  The  probe 
sampling  sequence  pertinent  to  this  experiment  lasts  a few  seconds.  It  is 
obvious  that  during  that  time  a steady-state  has  been  established  since  probe 
coverage  has  reached  one  molecular  layer  in  a matter  of  10  to  50  msec  in  our 
range  of  operation. 

Fournier  and  Pigache  have  reported  in  Ref.  57  a cleaning  sequence 
of  10  minutes  for  an  argon  beam  on  copper  or  stainless  steel  surfaces,  with 
the  following  properties: 


0.01  Coulomb/cm2  @ 2 KeV  ( 5 x 10' 


torr) 


This  represents  roughly  the  same  power  consumed  as  that  of  resistive  (ohmic) 
heating  at  500°C  for  the  same  period.  For  tungsten  and  for  the  region  of 
operation  typical  for  our  experimental  conditions,  we  can  expect  the  following 
estimates  (Refs.  79>8o) 


cleaning  time:  10  minutes 

beam  energy:  1 KeV  argon+ 

beam  density:  80  pamp/cm2 

We  have  taken  the  following  data  from  the  literature: 

(1)  Binding  energy  on  tungsten  (300°K) 

H2O , N25  O2 , CO2  — 100  Kcal/mole 

CO,  CH4,  K2  - 50  Kcal/mole 

Ar  m2  Kcal/mole 

(2)  Sticking  coefficient  on  tungsten  wire  (300°K) 

N2,  H2j  Ar  m 0.1  - 0.3 

02  m 0.2 

The  necessary  cleaning  power  reported  by  F-P  amounts  to 

_5 

Icn  cleaning:  8 x 10  KW/cmc  @ 10  min. 

Ohmic  cleaning:  10  x 10  5 KW/cm2  @ 10  min. 

In  fact,  we  can  see  that  some  power  is  being  delivered  to  the  surface  layers 
and  probably  suffices  to  knock  off  most  of  the  top  layers.  For  room  temperature 
in  fact,  at  5 x 10  5 torr,  the  N2  contributes  a non-negligible  iirpingement 
power  of  the  order  of 


We  estimate  the  number  of  foreign  molecules  stuck  on  the  surface 
of  the  probe  with  a sticking  coefficient  a = 0.2  on  tungsten  to  get  a rate 
at  a pressure  of  5 x 10-5  torr  (300°K)  of 

rsurf  -a  ro  - 4 x 1Ql5  cm’2  sec_1 

If  the  maximum  binding  energy  is  set  at  about  100  Kcal/mole,  this  yields  for 
each  molecule 


hind.  - 7 x 10'19  joule 


If  we  assume  that  the  ions  share  half  of  their  incoming  energy  and 
that  the  neutral  'hit'  on  the  surface  shares  also  half  with  its  neighbours, 
then  an  ion  delivers  roughly 


E E. 

clean  4 ion 

at  1 KeV  therefore,  one  requires  about 

T ~ 8 x 1013  ion/sec-cm2 
ion  ' 


or  a corresponding  power 


W.  ~ 1 x 10-5  KW/cm2 
ion  ' 

clean 

This  simplified  model  falls  a bit  short  of  the  figure  reported  by 
Pigache  for  copper  and  stainless  steel.  It  is  quite  plausible,  however,  to 
assume  that  the  incoming  ion  would  share  less  than  half  its  kinetic  energy 
and  that  moreover  a compound  sticking  coefficient  would  be  obtained  because 
of  the  many  species  accumulated  on  the  surface.  This  sinplified  model  is 
therefore  taken  as  accurate  to  within  an  order  of  magnitude.  The  impingement 
currents  required  (10-100  pA)  are  in  fact  easily  obtainable  with  the  standard 
operating  conditions  of  the  ion  source  used  for  this  experiment. 
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APPENDIX  D - MEASUREMENT  OF  ION  SPEED  RATIO 


USING  SMALL  COLLIMATED  FARADAY  CUPS 


/ 

A method  proposed  by  Fournier  and  Pigache  (Ref.  25)  outlines  a 
possibility  for  measuring  the  ion  speed  ratio  in  mesothermal  plasma  flows 
using  small  collimated  Faraday  cups.  The  basic  principles  invoked  are  not 
new.  However,  these  types  of  collectors,  used  mostly  as  ionospheric  plasma 
probes,  are  greatly  refined  by  achieving  an  improved  angular  sensitivity. 

A sinple  analytical  treatment  of  the  current  collection  mechanism  is  required 
to  obtain  the  basis  for  the  measuring  technique.  Briefly,  it  reduces  to 
using  the  fact  that  the  collimated  cup  is  made  velocity  selective  when  its 
attitude  with  respect  to  a plasma  stream  is  varied,  from  a parallel  to  a 
transverse  alignment  with  the  bulk  velocity  of  the  ions  in  the  flow.  The 
collimated  cups  have  the  distinct  advantage  over  the  cylindrical  probes  that 
the  current  collection  is  not  sensitive  to  the  presence  of  a randomized  slow- 
ion  content  within  the  beam.  This  comes  from  the  fact  that  the  'open  window' 
constitutes  a solid  angle  considerably  smaller  than  kir.  This  is  a situation 
virtually  the  opposite  of  the  one  encountered  with  the  'bare'  probe  which 
sees  nearly  kv  s'teradians. 

The  plasma  conditions  which  prevail  are  compatible  with  the  conditions 
pertinent  to  our  investigation.  The  electrons  are  randomized  and  possess  a 
given  temperature , smaller  than  say  a few  volts,  and  the  sheath  thickness  is 
large.  The  probe  orifice  is  considerably  smaller  than  the  Debye  length.  Thus 
the  authors  assume  that  no  electrostatic  field  effects  will  be  influencing  the 
ions  near  the  entrance  of  the  cup.  They  also  assume  that  the  space  within 
the  collecting  enclosure  is  field-free  and  that  no  significant  effects  are 
introduced  because  of  contamination  of  the  internal  surfaces  of  the  collector. 

The  system  is  shown  in  Fig.  A-2.  The  collector  is  left  floating,  for  simplicity, 
while  an  electron  suppressing  grid  is  biased  at  some  negative  voltage  - 9 volts 
in  this  case.  Therefore,  the  problem  is  reduced  to  a siirplified  'kinetic  theory' 
collection  of  the  ions  in  the  stream,  i 

The  cup  is  first  aligned  with  the  beam  geometric  centre-line  and  is 
pitched  about  its  first  aperture  by  an  angle  a,  from  0°  to  90°.  The  collector 
has  an  angular  acceptance  of  2(L(3) . The  ions  in  the  stream  have  a predominant 
directed  velocity  v and  a random  velocity  w.  The  system  of  coordinates  (spherical 
symmetry  w,  9,0)  is  chosen  such  that  the  x-axis  will  be  arbitrarily  aligned  with 
the  direction  of  v.  The  ion  distribution  away  from  the  cup  is  given  by  the 
equilibrium  Maxwellian  formulation,  as  was  assumed  throughout  this  investigation, 
i .e., 


f oc  exp(-[(w  - v)2  + w 2 + w 2]) 
w r *•  L v x ' y z -1 

The  directed  speed  is  given  by 

|v  | = u cosQi 


where  u is  the  average  bulk  ion  velocity  and  (a)  is  the  angle  between  the  geometric 
centreline  of  the  beam  and  the  x-axis.  Therefore,  the  current  density  is  given  by 
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4,  *e  J rw±d* 

open 

where  u is  the  velocity  along  some  line-of-sight  from  the  collector,  which  is 
at  an  angle  (£)  from  the  x-axis,  i.e. 

|u|  = |w|  sin0  cos(9  - P) 


The  integration  must  be  carried  out  over  the  entire  'OPEN  AREA*  viewed  by  the 
collector,  and  which  is  defined  by  the  double-aperture  system.  This  range  of 
integration  is  defined  using  the  following  limits: 

0 < w < oo 

p-Ap<e<p  + Ap 

I - A0  < 0 < |+ 


The  ion  current  to  the  collector  is  found  by  integrating  (j)  over  all  possible 
lines-of-sight  (3)  allowed  by  the  collimated  orifices.  Thus 


30, a) 


■f 

open 


dw  d0  d0  w3cos(8  - P)sin20  exp(-[w2  - 2wv  cos0  cos0  + v2]} 


where  v = v(a)  and 


I’ (a)  = JO,a)de 

v- 


A nondimensional  expression  is  obtained  by  normalizing  this  result  using  the 
value  at  zero  degree  of  incidence  (a  = 0) , i.e.  the  nondimensional  ion  current 
is  given  by 

^ ■ i^oy 

The  analytical  expression  for  j(p,a)  was  obtained  from  Ref.  25  and  a 
plot  of  1(a)  was  generated  by  numerical  integration  for  various  ion  speed 
ratios.  These  results  are  presented  in  Fig.  A-3.  By  cross-plotting  the  peak 
half-width  and  the  corresponding  speed  ratio  a master  plot  can  be  obtained  such 
as  the  one  shown  in  Fig.  A-4.  Such  a plot  readily  yields  the  relation  between 
01  and  Sj.  Consequently  the  range  of  sensitivity  of  the  method  can  be  assessed 
for  each  aperture  angle.  A measurement  of  the  absolute  peak  half- width  taken 
from  a current  versus  angle  trace  can  directly  yield  the  ion  speed  ratio.  Two 
experimental  cases  reproduced  in  Fig.  A-5  are  shown  for  comparison  in  Fig.  A-4. 
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The  speed  ratios  assumed  for  these  two  examples  are  about  44  and  34,  respectively, 
according  to  the  simple  geometrical  (conical  beam)  model  outlined  in  Chapter  III. 
The  measurements  are  short  of  these  model  values  by  50%  for  the  smaller  values  of 
Si  and  by  20 % for  the  large  Sj_.  Although  preliminary  results  show  a limited 
agreement  we  can  recognize  the  relatively  promising  features  of  this  measuring 
method.  The  relatively  simple  construction  has  the  rewarding  advantage  of 
giving  an  excellent  resolution  of  the  beam  density  which  can  prove  useful  for 
accurage  mapping  of  the  stream  (see  Fig.  III-6) . As  with  the  slender  cylindrical 
probe  used  in  this  study,  the  current  angular  response  of  the  cups  is  exceptionally 
sharp  for  small  angles  about  the  main  stream  direction.  This  makes  them  quite 
useful  in  locating  the  local  velocity  vector  within  the  beam.  Although  the 
currents  collected  by  the  cups  are  considerably  smaller  than  the  ones  achieved 
using  our  cylindrical  probes,  under  similar  operating  flow  conditions,  the 
signals  monitored  were  significantly  'cleaner',  as  we  have  mentioned  earlier, 
because  of  the  absence  of  a background  noise  (DC  offset)  coming  from  slow 
resident  ions  in  the  beam. 
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APPENDIX  E - CRITERIA  FOR  COLLISIONLESS  REGIME 


The  various  criteria  delineating  the  collisionless  plasma  regime, 
which  is  largely  satisfied  in  this  investigation  involving  the  use  of  cylindrical 
electrostatic  probes  in  a high  speed,  low  density  plasma  flow,  are  outlined 
schematically  below: 

(l)  The  plasma  is  said  to  be  rarefied  when  mechanisms  involving  close  encounters 
are  negligible  in  comparison  with  collective  phenomena,  i.e.  (macroscopic 
force- field  assumption) 


n A®  » 1 

(2)  In  the  case  where  condition  (l)  is  finite,  an  extra  condition  must  be 

realized  so  that  large  disturbances  caused  by  individual  particles  can  be 
neglected  and  it  reads 


Vp>>d 

where  d in  this  case  refers  to  the  probe  diameter. 

(3)  One  can  neglect  the  effects  of  neutrals  on  the  charged  particle  collection 
if 


\ » d,  X » d,  \ » d 

cl  Q 1 

provided  all  mean-free  paths  (mfp)  have  been  carefully  analyzed  (see 
Appendix  F) . 

(4)  Finally,  with  all  three  conditions  satisfied,  the  collector  must  be  a 

'probe'  in  the  true  sense,  such  that  the  sampling  currents  are  negligible 
coirpared  with  the  currents  responsible  for  maintaining  the  plasma  fluid 
(Ref.  13),  viz 


i . « 

probe 


I 


source 


and  that  the  probe  serves  as  an  ideal  sink  for  all  charged  species. 


APPENDIX  F - LENGTH  AND  TIME  SCALES 


(l)  Pertinent  Mean-Free-Paths 

The  various  mean-free-paths  (A)  can  be  described  schematically  as 
follows  for  a singly  ionized  argon  plasma,  where  e,  i,  a,  represent  electrons, 
ions  and  atoms,  respectively  and  P is  the  neutral  pressure 

(a)  X ~ 4.l8  P-1  [cm]  for  P [p.Hg] 

ClrCl 

Since  only  moderate  temperatures  are  involved  and,  if  we  assume  that  no 
significant  rate-process  occurs  ( charge- exchange , recombination. ..) , then  for 
the  plasma  regime  pertinent  to  this  study  we  have* 


/Mu.  c. 

A.  . ~ ( — > ■ - 

1a  = V Z,e 


~ A . 
ai 


and  for  charge-transfer , for  an  ion  energy  smaller  than  ~ 1.5  KeV 


Aia  — 3*3  x n”1  x 1014  [cm] 
For  an  electron  temperature  of  1 eV  then 


A ~ 1.2  n"1  x 1016  [cm] 

6a<  cl 


(d)  A a 2.5[n  in  A ]-1  • T 2 x 1013  [cm] 
ei  e eeJ  L J 


for  Tg  in  [ eV]  and  where 


A s 
ee 


■®[rVv 


Moreover  if  we  assume  that 


we  get 


v • - T,  V. 

ei  V M / le 


A.  2:  A . ( — = 

le  ei  V ad 


- — A . 

Oft'  61 


(g)  A.  ~ A 

' 7 ii  1 ee 


: ion  mobility . 


T.  \2  r in  A 


T.  r Xi. 

r)  • ;i 


Z4  in  A. 


(2)  Pertinent  Time  Scales 


Using  the  following  conditions  reported  herein, 


r “1;  Te-leV 
e J 


n.  < 108  cm"3 


n £ 5 x-1011  cm"3  (test  section) 


and  referring  to  (l)  we  can  .see  that  all  the  pertinent  mean-free-paths  (mfp) 
are  larger  than  the  te'st  section  diameter  (2.3m).  If  we  now  denote  the  character- 
istic  time-scales  by  t.  we  can  see  that  these  time- scales  are  all  significantly 
longer  than  the  ion  drift-time  past  the  probe  (tv^^^).  We  summarize  them  below. 


* t . 
ei 


\i  -5 

<7">  2:  5 x 10  sec. 
e 


* T. 

ie 


~ (M/m)  Tgi  ~ l.h  x 10'2  sec. 


~ (Pv^)"1  Tei  ~ 2 x 10-5  sec. 


- (M/m)2  t ~ 3 x 10” 4 sec. 
66 


Taa  ~ 10  sec' 


~ 1.8  x 10-4  sec. 


. CJ'TT  C 

rD  = ion  disturbance  time  = ~ 5 x 10  sec. 

Ti 


* Tcprift  = ion  drift- time 


(characteristic  length  £)  u 


_ l y -6 
= — <10  sec. 


Tmfr  = monolayer  formation  time  ~ 0.2  sec. 
(@  air,  10  torr) 


Ttrace  = characfceristic  'sample’  time  ~ few  seconds. 
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APPENDIX  G - TANGENT  METHODS  TO  MEASURE  Te 


In  order  to  simplify  the  measurements  of  the  electron  ten$?erature 
from  ideal  Langmuir  probe  characteristics,  a geometrical  construction  allows 
a probe  user  to  extract  an  approximated  value  directly  from  the  recorder  I-V 
trace,  for  both  single  and  double-probe  configurations.  Using  Sonin's  analysis 
(Ref.  13),  the  single-probe  collection  can  be  represented  by  the  relation  of 
the  net  current 


i = 3e  - ^i 


where 

j = jQexp(-ip)  (Boltzmann) 

and  ^ = e0/k T . Therefore  since  is  independent  of  Te 

( “e  V1  , _ 

di£ 


(^n  je) 


Thus  at  the  floating  potential  we  have  j 3 0 and 


kT 


1 dJil 

v-»  A 

1 |( 

Therefore,  using  the  simple  construction  shown  in  Fig.  A-6  for  the  single-probe 
characteristics,  we  have 


kT 


c + b 


= a 


such  that  a quick  estimate  of  Te  can  be  obtained  by  taking  the  tangent  at  the 
floating  potential  point. 


For  a double-probe  constructed  with  two  identical  probes  we  can  use 
the  ideal  probe  characteristic  reproduced  in  Fig.  A-6.  If  (l)  and  (2)  indicate 
the  collection  regions  for  probe  No.  1 and  probe  No.  2,  respectively,  then 
from  Ref.  66  we  have 


Ja.  = J* 


exp 


e 


kT 


- 1 


= J (1) 


= J* 


1 - exp 


30] 


= J (2) 


G-l 


where  j*  j.s  the  constant  i 


on  current  to  the  probes  obtained  by  extrapolation 

J 


- e 0 

J*  = c exp  I space 


Thus  if  + At>s  = o (ideally) 


kT 


e _ 


2 


t$l 


S.iLS'S-  *»  «-  ™uase 

0 - o,  and  the  extrapolated  values  of  j*,  i?e.  * slope’  °r  t£mgent  at 


APPENDIX  H - PROBE  STRUCTURAL  EFFECTS 


The  long  cylindrical  probes  are  treated  here  as  cantilevered  cylindrical 
thin  beams.  Indeed,  a number  of  structural  effects  could  threaten  their  integrity 
as  ideal  collectors.  Such  effects  can  be  summarized  as  follows: 

. flexure  or  bending  from  the  action  of  gravity, 

. oscillation  at  the  natural  frequency  or  at  induced  frequencies 
from  different  sources  such  as  probe  traversing  mechanism  and 
foundation  noise, 

. thermal  expansion, 

. non- uniformity  of  cross-section  caused  in  part  by  inherent 

metal  properties  or  by  deformation  occurring  in  construction: 
roughness,  extrusions,  kinks. 

After  inspection,  we  find  that  all  but  the  first  can  be  considered  negligible. 

The  pure  bending  relationships  are  given  by  (Ref.  83): 


(1) 


. _ Wi3  . . „ ^ ^ 

W = BeT  = rca-x-  tip  deflection 


where  E is  Young's  Modulus  and  W is  the  weight 


_ _ rt£ 


W = p 


T 


and  the  moment  of  inertia  is  given  by 


I e 


_ 7rd 


~r 


For  small  angles  then, 


_ 4 


(2)  0s  ~ 'i  l j _ initial  tip  slope  @ max,  deflection 


(3) 


= 3-515 


nat 


El 

gj?  ( — ) s natural  frequency 


Upon  rearranging  (l)  we  get  the  normalized  maximum  deflection 


max 


^max  ) _ 4piR3 

d ' E 


where  R is  the  probe  aspect  ratio  (i/d)  and  (p)  is  the  material  density  and 


H-l 


9 = artan 

s 


i I* 

4 max 

3 R 


) 


Wnat  ~ 3,515^^7§0Hd  * (*5ax) 

for  (d)  in  mils  and  g = 32.2  ft/sec2.  Using  a typical  range  of  parameters 
pertinent  to  this  study,  we  see  that  one  obtains  significant  bending  for  a 
long  slender  tungsten  probe,  e.g.,  using 


R = 500  and  d = 0.004  in. 


then 


I* 

max 


~ 6.5 


0 

s 


~ 1° 


u 


nat 


~ 30  Hz 


Therefore  some  care  must  be  taken  with  our  long  slender  probe  (R  ~ 472.5, 
d = 0.004  in.)  since  anomalous  probe  effects  can  be  perceived,  along  with 
oscillations,  which  are  within  the  sensitivity  of  the  equipment.  In  fact, 
visual  observation  indicated  that  such  oscillations  were  present  when  a 
nearby  experimental  facility  was  operated  (UTIAS  Low  Density  Wind  Tunnel), 
and  measurements  were  restricted  to  'quiet'  periods  when  using  the  longer 
probes. 
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APPENDIX  I - INVERSE  PROBLEM  FOR  LONG  CYLINDRICAL  PROBES 


A systematic  method  outlined  by  Sajben  and  Lee  (Ref.  84)  demonstrates 
the  concept  of  solving  the  inverse  problem  for  the  collection  of  electrostatic 
probes  in  general,  i.e.  when  the  analytical  description  of  the  current  collection 
is  known  a priori,  one  can  uniquely  determine  the  fundamental  state  properties 
of  the  plasma.  The  interesting  advantage  of  the  method  consists  in  the  deter- 
mination of  the  plasma  properties  ne,  Te,  0Space  f°r  each  point  on  the  current- 
voltage  characteristic  such  that  a single  trace  can  yield  a set  of  statistically 
redundant  measurements  of  these  quantities. 


In  order  to  verify  the  overall  ability  of  bhe  method  we  have  attempted 
to  apply  it  to  the  steady- dynamic  case  using  the  analytical  description  obtained 
in  Chapter  I,  together  with  some  experimental  I-V  traces.  Although  no  quantitative 
description  is  provided  here  the  preliminary  results  have  shown  that  the  scheme 


yields  two  solutions  for  Te  and  that  the  region  of  sensitivity  is  reduced  to  a 


very  small  range  of  potentials  around  the  floating  value  (Vp  ± 5 volt) . 


This  fact  that  the  solution  for  Te  is  not  unique  and  that  the  sensitiv- 
ity of  the  method  is  quite  limited  is  interesting,  as  it  corroborates  both  the 
relative  non-sensitivity  of  the  ion  current  collection  by  the  parallel  probe  to 
the  electron  temperature,  and  the  duality  of  the  solution  also  obtained  by  our 
approximation  described  in  Appendix  J.  Therefore,  we  conclude  that  this  method 
would  be  of  no  significant  help  for  treating  the  dynamic  case  involving  long 
negative  parallel  cylinders,  although  it  can  prove  useful  to  probe  users  for  the 
study  of  quiescent  plasmas. 


APPENDIX  J - STATUS  OF  THE  MASS-MOTION  EFFECT 


The  study  of  the  end  effect  has  actually  contributed  to  shed  more 
light  on  the  long  debated  question  whether  or  not  the  bulk  displacement  of  the 
plasma,  or  the  ’mass-motion',  would  effect  the  characteristics  of  the  current 
collection  of  cylindrical  Langmuir  probes.  We  provide  in  Fig.  A-7  an  evaluation 
of  the  current  status  on  this  subject  in  a fashion  advocated  by  French  (Ref.  7). 

Two  flow  regimes  are  considered,  namely  that  of  collisionless  thin  and 
thick  sheath,  respectively* 

(1)  Kn  » 1;  A » d » 

(2)  Kn  » 1;  A » A^  > d 

The  second  case  corresponds  to  the  regime  of  interest  to  this  study,  and  we  have 
indicated  it  by  the  hatched  region  in  Fig.  1-1.  We  indicate  by  YES  that  either 
a theory  or  an  experimental  investigation,  or  both,  confirm  the  existence  of  an 
effect  due  to  the  bulk  displacement  of  the  plasma  relative  to  the  probe.  A NO 
indicates  that  no  mass-motion  effect  is  observed  or  predicted.  An  * indicates 
that,  to  our  knowledge,  the  information  pertinent  to  the  corresponding  conditions 
is  incomplete  in  the  literature. 

The  region  implemented  by  the  results  obtained  in  this  investigation, 
using  the  analysis  of  J.  R.  Sanmartin,  is  enhanced  by  the  heavy- line  frame  for 
the  ion  collecting  region.  Another  area  of  interest  was  covered  recently  by 
the  work  of  Taillet  et  al  (Ref.  2k).  They  actually  demonstrated  both  analytically 
and  by  conputer  simulation  that  the  electron  collection  region  would  indeed  be 
altered  by  the  plasma  motion.  This  interesting  prediction  has  effectively  been 
corroborated  by  the  experimental  results  obtained  by  Koopman  and  Segall  using  a 
high-speed  plasma  produced  when  a pulsed  laser  beam  impacts  on  copper  or  aluminum 
targets  (Ref.  32). 

Finally,  it  can  be  pointed  out  that  the  experimental  apparatus  used  in 
conjunction  with  the  present  investigation  could  prove  useful  in  conpleting  this 
survey  for  the  regimes  requiring  more  elucidation. 


*Kn: Knudsen  number  = A/r  and  r the  probe  radius  (d  = 2r  ) 
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APPENDIX  K - ESTIMATE  OF  ELECTRON  TEMPERATURE 


As  we  have  intuitively  indicated  in  Chapter  I the  end  effect  is 
relatively  insensitive  to  the  elec£ron  temperature.  However,  we  have  seen 
that  the  effective  scaling  length  6 of  the  potential  field  around  the  probe 
needs  to  be  adjusted  because  the  finite  probe  draws  ions  from  a region  within 
the  cylindrical  sheath  which  is  smaller  than  the  one  estimated  by  the  classical 
sheath  theory.  This  restricted  region  is  effectively  defined  by  the  limits: 


Ti 


< 3 


This  has  justified  the  definition  of  a characteristic  field  length  6 which  is 
described  by  Sanmartin  using  the  following  exact  relation 


6"1  s <2n(e_1)  + Y(e,^p)  (l) 

Sanmartin,  in  fact,  does  not  define  Y(e,ip  ) analytically,  but  rather  provides 
the  graphical  description  we  have  reproduced  in  Fig.  A-8.  We  have  fitted 
these  curves  numeri cally . This  has  yielded  the  following  estimate  (±  %) 

Y(e,^p)  ~ 0.23  - o.l  ^n(e-1)  + 0.42|^p|°’3  (2) 

This_estimate  was  compared  with  the  approximation  explicitly  given  by  Sanmartin 
for  6 in  1.39*  Therefore,  using  Eq.  1.38  his  results  read 


6”1  m in(e-1)  - 0.5  + i ny  2^p6  (3) 

whdreas 

^fit  ~ 0,9  -M6'1)  + 0.23  + 0.42  ^°“3  (4) 

This  comparison  is  outlined  in  two  tables  shown  in  Fig.  A- 9 for  a range  of 
parameters.  The  agreement  is  excellent  except  for  the  comparison  between  the 
values  of  in(e  x)  predicted  at  large  probe  bias  (ip  ). 

P 

Since  our  self-consistent  method  described  in  Section  IV. 2.3  actually 
determines  5 using  experimental  measurements  exclusively  related  to  ion  features 
and  to  probe  characteristics,  it  is  conceivable  to  expect  some  electron  temper- 
ature information  out  of  the  value  of  6sepf.  The  relatively  small  sensitivity 
of  6 to  electron  tenperature  is,  however,  apparent  from  Fig.  A- 11.  We  have 
plotted  the  expression  obtained  in  Eq.  (3)  using  the  exact  relation  for  e and 
^p,  namely 


K-l 


e = 


W 


- VV  = H6 

e 


and  we  notice  that  a three-order  of  magnitude  change  in  Te  implies  merely  a 
± 20%  change  in  5,  using  a typical  set  of  experimental  conditions. 

Nevertheless,  an  attempt  to  measure  Te  indirectly  was  made.  This 
method  implies  that  once  6 is  known 


6 = 5 


self 


one  must  find  an  agreement  between  the  value  of  function  Y(e,ip  ) obtained 
from  (l)  and  from  (2),  that  is:  ^ 

Yi  = Ya(e,^p)  = ®gglf  - Me"1) 

Y2  = Y z{e,ip  ) ~ 0.23  - 0.1  ^n(e-1)  + 0.42  t/P'3 
P P 

with 


F = |Yi  - Y2 | = 0 


(5) 


If,  however,  one  plots  F(e),  one  notices  that  for  certain  values  of  5 condition 
(5)  is  never  satisfied.  This  has  led  us  to  use  an  approximate  solution  of  (5) 
using 


£-o 


@ e = e 


m 


as  shown  in  Fig.  A- 10.  The  estimates  of  e using  this  scheme,  and  labelled  %, 
are  inherently  only  marginally  accurate.  For  instance,  an  error  of  less  than 
10%  on  the  value  of  in(e  1)  leads  to  error  of  about  50%  on  the  value  of  e with 
e ~ 0.01.  Nevertheless,  an  average  value  of  Te  was  obtained  with  this  relatively 
coarse  method.  It  is  on  average  about  100%  in  disagreement  with  the  direct 
measurements  of  Te  (Chapter  III),  i.e. 

Te*=®e 


as  shown  in  the  table  presented  in  Fig.  A-ll.  This  is  nevertheless  quite 
interesting  since  the  single-probe  results  also  suggest  a measured  value  of  Te 
about  a factor  of  2 too  high.  This,  in  fact,  constitutes  an  extra  indication 
that  the  bulk  motion  does  introduce  some  distortion  in  the  electron  temperature 
sensitive  region  of  the  current  versus  voltage  probe  characteristic,  as  was 
suggested  in  typical  I-V  probe  responses  (see  Fig.  IV-2) , and  as  was  corroborated 
elsewhere  (Appendix  j) . 
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